library(lmerTest)
library(ggeffects)
library(dplyr)
library(report)
package ‘report’ was built under R version 4.0.5
library(r2glmm)
fullTest <- read.csv("../Cleaning/output/fullTest.csv")
fullTrain <- read.csv("../Cleaning/output/fullTrain.csv")

traitsFreqs <- read.csv("../Cleaning/output/traitFreqOverUnder.csv")
traitsFreqs <- rename(traitsFreqs, props = optionChoiceN)
fullTest <- merge(fullTest, traitsFreqs[c("trait","props")], by = "trait")
fullTest$propCorr <- ifelse(fullTest$Estimator=="Underestimator", 1-fullTest$props, fullTest$props)

uSubs <- unique(fullTest$subID)

indDiffs <- fullTest[!duplicated(fullTest$subID),]
fullTest$ingChoiceN <- as.factor(fullTest$ingChoiceN)
fullTest$novel <- as.factor(fullTest$novel)
fullTest$selfResp.Z <- scale(fullTest$selfResp)
fullTest$SE.Z <- scale(fullTest$SE)
fullTest$iSE.Z <- scale(fullTest$iSE)
fullTest$oSE.Z <- scale(fullTest$oSE)
fullTest$predicted.Z <- scale(fullTest$predicted)
fullTest$slope.Z <- scale(fullTest$slope)
fullTest$entropy.Z <- scale(fullTest$entropy)
fullTest$WSR.Z <- scale(fullTest$WSR)
fullTest$neighAveOutSE.Z <- scale(fullTest$neighAveOutSE)
fullTest$neighAveAllSE.Z <- scale(fullTest$neighAveAllSE)
fullTest$neighAveInSE.Z <- scale(fullTest$neighAveInSE)
fullTest$novel <- as.factor(fullTest$novel)
levels(fullTest$novel) <- list("Trained"  = "0", "Held Out" = "1")
# prop.test(traitsFreqs$optionChoiceN, traitsFreqs$N, p=rep(.5,length(traitsFreqs$N)))
# 
# m <- glmer( ingChoiceN ~ trait + ( 1 | subID) , data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
#                                     optCtrl = list(maxfun = 100000)),
#     nAGQ = 1)
# 
# fullTest$trait <- as.factor(fullTest$trait)
# contrasts(fullTest$trait) <- contr.sum(148)
# m <- glm(optionChoiceN ~ trait, family = binomial,
#           data = fullTest
#          )
# summary(m)
summary(m)

Call:
glm(formula = optionChoiceN ~ 1, family = binomial, data = traitsFreqs)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
-0.81807  -0.21571   0.02476   0.21563   0.82291  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)
(Intercept)   0.1913     0.1652   1.158    0.247

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 13.212  on 147  degrees of freedom
Residual deviance: 13.212  on 147  degrees of freedom
AIC: 203.17

Number of Fisher Scoring iterations: 3
propMatrix <- matrix(nrow=148,ncol=7)
for(i in 1:148){
    traitDf <- subset(fullTest, Idx==i)
    test <- t.test(as.numeric(traitDf$ingChoiceN)-1, mu=.50)
    propMatrix[i, ] <- c(i, test$statistic, test$p.value, test$conf.int, test$estimate, test$parameter)
}
colnames(propMatrix) <- c("Idx", "stat", "p", "LCI", "UCI", "est", "param")
propMatrix <- as.data.frame(propMatrix)
propMatrix$trait <- traitsFreqs$trait
propMatrix <- propMatrix[order(propMatrix$p),]
library(corrr)
x <- indDiffs %>% 
    select(groupHomoph, DS:SING.Inter) %>%
  correlate() %>% 
  focus(groupHomoph) %>%
    arrange(groupHomoph)

Correlation method: 'pearson'
Missing treated using: 'pairwise.complete.obs'
x %>% 
  mutate(rowname = factor(rowname, levels = rowname[order(groupHomoph)])) %>%  # Order by correlation strength
  ggplot(aes(x = rowname, y = groupHomoph)) +
    geom_bar(stat = "identity") +
    ylab("Correlation Coefficient") +
    xlab("Individual Differences") + theme_grey(base_size = 9)  + theme(axis.text.x = element_text(angle = 90,hjust = 1)) +
  theme(axis.text.x = element_text( 
                           size = 9, angle = 45, vjust = 1)) + theme(axis.title.x = element_text(vjust=1.9)) + theme(axis.text=element_text(size=9),
        axis.title=element_text(size=9,face="bold")) + theme(legend.text = element_text(size=9)) + theme(panel.border = element_rect(colour = "black", fill = NA, size =1)) + theme(legend.title = element_blank()) + theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),
panel.background = element_blank(), axis.line = element_line(colour = "black"))

Familiarity predicts Reaction Time

summary(m)
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: log(RT) ~ fam + propCorr + desirability + (fam | subID) + (1 |      trait)
   Data: fullTest

REML criterion at convergence: 11321

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-8.3287 -0.5140 -0.0934  0.4175  7.9841 

Random effects:
 Groups   Name        Variance  Std.Dev. Corr
 trait    (Intercept) 8.900e-03 0.094340     
 subID    (Intercept) 4.609e-01 0.678881     
          fam         2.238e-05 0.004731 1.00
 Residual             7.821e-01 0.884369     
Number of obs: 4292, groups:  trait, 148; subID, 29

Fixed effects:
               Estimate Std. Error         df t value Pr(>|t|)  
(Intercept)   1.742e-01  2.418e-01  1.360e+02   0.720   0.4727  
fam           6.877e-03  3.544e-03  1.169e+02   1.940   0.0548 .
propCorr      8.663e-02  8.951e-02  4.223e+03   0.968   0.3332  
desirability -2.723e-02  3.800e-02  1.468e+02  -0.717   0.4747  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) fam    prpCrr
fam          0.172              
propCorr    -0.175  0.030       
desirabilty -0.776 -0.393 -0.023
optimizer (nloptwrap) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')

Replication of prior self-anchoring findings: Self-evaluations predicting ingroup evaluations

summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ selfResp.Z + propCorr + desirability + (selfResp.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  3050.2   3097.0  -1517.1   3034.2     2548 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.2611 -0.8119  0.4028  0.7411  4.3393 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.00000  0.0000       
 subID  (Intercept) 0.47448  0.6888       
        selfResp.Z  0.09608  0.3100   0.61
Number of obs: 2556, groups:  trait, 148; subID, 29

Fixed effects:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept)  -1.45628    0.62618  -2.326   0.0200 *  
selfResp.Z    0.27625    0.07804   3.540   0.0004 ***
propCorr      5.00392    0.32287  15.498   <2e-16 ***
desirability -0.12952    0.10317  -1.255   0.2093    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slfR.Z prpCrr
selfResp.Z   0.183              
propCorr    -0.204  0.052       
desirabilty -0.942 -0.106 -0.054
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')

Does similarity-weighted self-evaluation average predict ingroup choices?

summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ WSR.Z + propCorr + desirability + (WSR.Z | subID) +      (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  5030.1   5080.9  -2507.1   5014.1     4213 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-5.6590 -0.7921  0.3796  0.7276  3.6805 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000    0.000        
 subID  (Intercept) 5.967    2.443        
        WSR.Z       6.912    2.629    0.40
Number of obs: 4221, groups:  trait, 148; subID, 29

Fixed effects:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept)  -1.99361    0.69564  -2.866  0.00416 ** 
WSR.Z         0.74295    0.54509   1.363  0.17288    
propCorr      4.83641    0.26082  18.543  < 2e-16 ***
desirability -0.11835    0.08155  -1.451  0.14672    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) WSR.Z  prpCrr
WSR.Z        0.289              
propCorr    -0.158  0.020       
desirabilty -0.669 -0.087 -0.053
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')

Does self-evaluation weighted similarity predict ingroup choices?

Do cross-validated similarity*self-evaluation predictions predict ingroup choices?

Do cross-validated similarity*self-evaluation predictions predict ingroup choices, regardless of whether it was seen prior or not?

m <- glmer( ingChoiceN ~ predicted.Z * novel + ( predicted.Z + novel | subID) , data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ predicted.Z * novel + (predicted.Z + novel | subID)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4836.3   4898.4  -2408.2   4816.3     3680 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.2577 -1.0435  0.6154  0.8642  2.1657 

Random effects:
 Groups Name          Variance Std.Dev. Corr       
 subID  (Intercept)   0.297398 0.54534             
        predicted.Z   0.015375 0.12400  -0.95      
        novelHeld Out 0.004533 0.06733   0.10  0.22
Number of obs: 3690, groups:  subID, 29

Fixed effects:
                           Estimate Std. Error z value Pr(>|z|)   
(Intercept)                0.309950   0.116300   2.665   0.0077 **
predicted.Z                0.259248   0.107170   2.419   0.0156 * 
novelHeld Out              0.002235   0.075330   0.030   0.9763   
predicted.Z:novelHeld Out -0.053348   0.077558  -0.688   0.4915   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z nvlHlO
predicted.Z -0.321              
novelHeldOt -0.308  0.087       
prdctd.Z:HO  0.067 -0.289 -0.043
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
tidy(m,conf.int=TRUE,exponentiate=TRUE,effects="fixed")
r2beta(m)
boundary (singular) fit: see help('isSingular')
boundary (singular) fit: see help('isSingular')
boundary (singular) fit: see help('isSingular')
boundary (singular) fit: see help('isSingular')
                     Effect   Rsq upper.CL lower.CL
1                     Model 0.012    0.021    0.007
2               predicted.Z 0.008    0.015    0.004
4 predicted.Z:novelHeld Out 0.000    0.002    0.000
3             novelHeld Out 0.000    0.001    0.000
ggpredict(m, c("predicted.Z", "novel")) %>% plot(show.title=F)+ xlab("Cross-Validated Self-Descriptiveness") + ylab("Likelihood of Ingroup Choice") + jtools::theme_apa()  + scale_color_discrete(labels = c("Trained","Held-Out"))
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.
Scale for 'colour' is already present. Adding another scale for 'colour', which will replace the existing scale.
ggsave("~/Documents/UC Riverside/Studies/Self-Anchoring/Figures/CVpredictionGeneralize.tiff",dpi=600)
Saving 7.29 x 4.51 in image

Does generalization depend on outdegree?

m <- glmer( ingChoiceN ~ predicted.Z * novel * outDegree + ( predicted.Z + novel | subID) , data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ predicted.Z * novel * outDegree + (predicted.Z +      novel | subID)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4832.9   4919.9  -2402.5   4804.9     3676 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.3855 -1.0379  0.5867  0.8615  2.2460 

Random effects:
 Groups Name          Variance Std.Dev. Corr       
 subID  (Intercept)   0.29679  0.54479             
        predicted.Z   0.01238  0.11127  -0.91      
        novelHeld Out 0.00446  0.06678   0.09  0.33
Number of obs: 3690, groups:  subID, 29

Fixed effects:
                                     Estimate Std. Error z value Pr(>|z|)  
(Intercept)                          0.243639   0.138851   1.755   0.0793 .
predicted.Z                          0.291423   0.134287   2.170   0.0300 *
novelHeld Out                       -0.148255   0.142113  -1.043   0.2968  
outDegree                            0.003400   0.003735   0.910   0.3627  
predicted.Z:novelHeld Out           -0.331784   0.143766  -2.308   0.0210 *
predicted.Z:outDegree               -0.002504   0.003750  -0.668   0.5043  
novelHeld Out:outDegree              0.006894   0.005731   1.203   0.2290  
predicted.Z:novelHeld Out:outDegree  0.013065   0.005807   2.250   0.0245 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z nvlHlO outDgr pr.Z:HO pr.Z:D nvHO:D
predicted.Z -0.184                                           
novelHeldOt -0.442  0.064                                    
outDegree   -0.555  0.010  0.544                             
prdctd.Z:HO  0.058 -0.430  0.026 -0.056                      
prdctd.Z:tD  0.058 -0.577 -0.054 -0.125  0.538               
nvlHldOt:tD  0.364 -0.040 -0.847 -0.649 -0.041   0.082       
prdc.Z:HO:D -0.035  0.371 -0.038  0.080 -0.844  -0.646  0.040
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
tidy(m,conf.int=TRUE,exponentiate=TRUE,effects="fixed")
r2beta(m)
boundary (singular) fit: see help('isSingular')
boundary (singular) fit: see help('isSingular')
boundary (singular) fit: see help('isSingular')
boundary (singular) fit: see help('isSingular')
                               Effect   Rsq upper.CL lower.CL
1                               Model 0.013    0.023    0.008
2                         predicted.Z 0.003    0.007    0.000
5           predicted.Z:novelHeld Out 0.002    0.005    0.000
8 predicted.Z:novelHeld Out:outDegree 0.001    0.005    0.000
7             novelHeld Out:outDegree 0.000    0.003    0.000
3                       novelHeld Out 0.000    0.002    0.000
4                           outDegree 0.000    0.002    0.000
6               predicted.Z:outDegree 0.000    0.002    0.000
ggpredict(m, c("predicted.Z", "outDegree" ,"novel")) %>% plot(show.title=F)+ xlab("Cross-Validated Self-Descriptiveness") + ylab("Likelihood of Ingroup Choice") + jtools::theme_apa()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.
ggsave("~/Documents/UC Riverside/Studies/Self-Anchoring/Figures/CVpredictionGeneralize.tiff",dpi=600)
Saving 7.29 x 4.51 in image

Does generalization depend on indegree?

Nope

summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ predicted.Z * novel * inDegree + (predicted.Z +      novel | subID)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4842.6   4929.6  -2407.3   4814.6     3676 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.3708 -1.0425  0.6081  0.8478  2.1489 

Random effects:
 Groups Name          Variance Std.Dev. Corr       
 subID  (Intercept)   0.297726 0.54564             
        predicted.Z   0.014331 0.11971  -0.94      
        novelHeld Out 0.004817 0.06941   0.06  0.29
Number of obs: 3690, groups:  subID, 29

Fixed effects:
                                     Estimate Std. Error z value Pr(>|z|)  
(Intercept)                         0.3072752  0.1489679   2.063   0.0391 *
predicted.Z                         0.2147101  0.1427624   1.504   0.1326  
novelHeld Out                      -0.0106251  0.1574787  -0.067   0.9462  
inDegree                            0.0001678  0.0045112   0.037   0.9703  
predicted.Z:novelHeld Out          -0.1367979  0.1585589  -0.863   0.3883  
predicted.Z:inDegree                0.0019902  0.0043923   0.453   0.6505  
novelHeld Out:inDegree              0.0005728  0.0065070   0.088   0.9298  
predicted.Z:novelHeld Out:inDegree  0.0040134  0.0066092   0.607   0.5437  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z nvlHlO inDegr pr.Z:HO pr.Z:D nvHO:D
predicted.Z -0.172                                           
novelHeldOt -0.495  0.045                                    
inDegree    -0.626 -0.007  0.598                             
prdctd.Z:HO  0.043 -0.485  0.019 -0.033                      
prdctd.Z:nD  0.041 -0.652 -0.035 -0.069  0.585               
nvlHldOt:nD  0.437 -0.020 -0.878 -0.693 -0.027   0.048       
prdc.Z:HO:D -0.024  0.429 -0.025  0.046 -0.873  -0.666  0.020

Neighboring Dependencies Predicting Choices

Generalization of Outdegree Neighboring Self-Evaluations

summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ neighAveOutSE.Z * novel + desirability + propCorr +  
    (neighAveOutSE.Z + novel | subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4419.5   4500.0  -2196.7   4393.5     3601 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.1997 -0.8726  0.4282  0.7774  4.3536 

Random effects:
 Groups Name            Variance  Std.Dev.  Corr     
 trait  (Intercept)     8.505e-15 9.222e-08          
 subID  (Intercept)     3.667e-01 6.056e-01          
        neighAveOutSE.Z 3.430e-02 1.852e-01 0.44     
        novelHeld Out   2.598e-02 1.612e-01 0.13 0.95
Number of obs: 3614, groups:  trait, 147; subID, 25

Fixed effects:
                              Estimate Std. Error z value Pr(>|z|)    
(Intercept)                   -2.05530    0.51479  -3.992 6.54e-05 ***
neighAveOutSE.Z                0.23453    0.07600   3.086  0.00203 ** 
novelHeld Out                  0.02107    0.08187   0.257  0.79692    
desirability                  -0.01230    0.08333  -0.148  0.88266    
propCorr                       4.79310    0.26952  17.784  < 2e-16 ***
neighAveOutSE.Z:novelHeld Out  0.03691    0.08016   0.460  0.64520    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) nAOSE. nvlHlO dsrblt prpCrr
nghAvOtSE.Z  0.105                            
novelHeldOt -0.057  0.167                     
desirabilty -0.931 -0.071  0.015              
propCorr    -0.227  0.045  0.006 -0.035       
ngAOSE.Z:HO -0.003 -0.375  0.013  0.007  0.032
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')

Does entropy (i.e., uncertainty) predict likelihood of ingroup choices?

Does a linear trend of similarity-based probabilities predict ingroup choices?

Does a linear trend of similarity-based probabilities predict ingroup choices, controlling for self-descriptiveness?

m <- glmer( ingChoiceN ~ scale(slope) + selfResp.Z  + ( scale(slope) + selfResp.Z | subID) + (1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ scale(slope) + selfResp.Z + (scale(slope) + selfResp.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  3324.8   3383.3  -1652.4   3304.8     2546 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.0960 -1.0011  0.5507  0.8387  2.3694 

Random effects:
 Groups Name         Variance Std.Dev. Corr       
 trait  (Intercept)  0.00000  0.0000              
 subID  (Intercept)  0.38198  0.6180              
        scale(slope) 0.04406  0.2099   -0.99      
        selfResp.Z   0.13234  0.3638    0.55 -0.66
Number of obs: 2556, groups:  trait, 148; subID, 29

Fixed effects:
             Estimate Std. Error z value Pr(>|z|)  
(Intercept)   0.27038    0.12721   2.126   0.0335 *
scale(slope)  0.21730    0.11988   1.813   0.0699 .
selfResp.Z    0.19319    0.08783   2.200   0.0278 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(s)
scale(slop) -0.406       
selfResp.Z   0.403 -0.454
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
tidy(m,conf.int=TRUE,exponentiate=TRUE,effects="fixed")

Does a non-parametric trend of similarity-based probabilities predict ingroup choices, controlling for self-descriptiveness?

m <- glmer( ingChoiceN ~ scale(nlslope) + selfResp.Z +  ( scale(nlslope) + selfResp.Z | subID) + (1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ scale(nlslope) + selfResp.Z + (scale(nlslope) +      selfResp.Z | subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  3327.3   3385.7  -1653.6   3307.3     2546 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.0256 -0.9895  0.5363  0.8481  2.3016 

Random effects:
 Groups Name           Variance Std.Dev. Corr       
 trait  (Intercept)    0.0000   0.0000              
 subID  (Intercept)    0.4056   0.6369              
        scale(nlslope) 0.2145   0.4632   -0.24      
        selfResp.Z     0.1390   0.3728    0.62 -0.44
Number of obs: 2556, groups:  trait, 148; subID, 29

Fixed effects:
               Estimate Std. Error z value Pr(>|z|)  
(Intercept)     0.35811    0.14766   2.425   0.0153 *
scale(nlslope)  0.17096    0.14329   1.193   0.2328  
selfResp.Z      0.20564    0.08781   2.342   0.0192 *
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(n)
scal(nlslp) -0.073       
selfResp.Z   0.393 -0.389
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
tidy(m,conf.int=TRUE,exponentiate=TRUE,effects="fixed")

Does a linear trend of similarity-based probabilities predict ingroup choices?

m <- glmer( ingChoiceN ~ slope.Z  + predicted.Z + ( slope.Z + predicted.Z | subID) + (1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ slope.Z + predicted.Z + (slope.Z + predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4653.3   4715.5  -2316.7   4633.3     3680 

Scaled residuals: 
   Min     1Q Median     3Q    Max 
-6.507 -1.001  0.472  0.811  2.610 

Random effects:
 Groups Name        Variance Std.Dev. Corr       
 trait  (Intercept)  0.0000  0.0000              
 subID  (Intercept)  0.4779  0.6913              
        slope.Z     20.1334  4.4870   -0.59      
        predicted.Z 21.3954  4.6255    0.57 -1.00
Number of obs: 3690, groups:  trait, 148; subID, 29

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)   0.2051     0.1384   1.481   0.1385  
slope.Z       1.5056     0.9093   1.656   0.0978 .
predicted.Z  -1.3237     0.9392  -1.409   0.1587  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z
slope.Z     -0.573       
predicted.Z  0.537 -0.994
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
tidy(m,conf.int=TRUE,exponentiate=TRUE,effects="fixed")

Backwards solution: Can you predict self-evaluations from similarity to ingroup and outgroup choices?

summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ eSE + (eSE | subID)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4785.1   4816.1  -2387.5   4775.1     3635 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.6814 -1.0139  0.5896  0.8615  2.8232 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 subID  (Intercept) 24.66    4.965         
        eSE         16.76    4.094    -1.00
Number of obs: 3640, groups:  subID, 25

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)
(Intercept)  -0.4974     1.1746  -0.424    0.672
eSE           0.6004     0.9522   0.631    0.528

Correlation of Fixed Effects:
    (Intr)
eSE -0.996
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: ingChoiceN ~ sSE * novel + (sSE | subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4791.8   4841.4  -2387.9   4775.8     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.6486 -1.0125  0.5929  0.8627  2.7918 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept)  0.000   0.000         
 subID  (Intercept)  1.134   1.065         
        sSE         26.117   5.111    -0.90
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)
(Intercept)  0.03974    0.29415   0.135    0.893
sSE          0.89300    1.28620   0.694    0.487
novel1       0.12039    0.32510   0.370    0.711
sSE:novel1  -0.36572    1.26144  -0.290    0.772

Correlation of Fixed Effects:
           (Intr) sSE    novel1
sSE        -0.928              
novel1     -0.427  0.366       
sSE:novel1  0.413 -0.372 -0.976
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(oSE) + (scale(oSE) | subID) + (1 |      trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  5482.7   5520.8  -2735.3   5470.7     4215 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.6879 -0.9939  0.5179  0.8500  2.8947 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.0000   0.0000       
 subID  (Intercept) 0.3871   0.6222       
        scale(oSE)  0.1068   0.3267   0.59
Number of obs: 4221, groups:  trait, 148; subID, 29

Fixed effects:
            Estimate Std. Error z value Pr(>|z|)  
(Intercept)  0.30322    0.12022   2.522   0.0117 *
scale(oSE)   0.08435    0.06924   1.218   0.2232  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
           (Intr)
scale(oSE) 0.504 
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(RSE) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(RSE) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4449.7   4498.6  -2216.8   4433.7     3364 

Scaled residuals: 
   Min     1Q Median     3Q    Max 
-2.284 -1.052  0.623  0.872  2.044 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr 
 trait  (Intercept) 2.669e-16 1.634e-08      
 subID  (Intercept) 2.409e-01 4.908e-01      
        predicted.Z 7.502e-05 8.662e-03 -1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                       Estimate Std. Error z value Pr(>|z|)   
(Intercept)             0.24110    0.11447   2.106  0.03518 * 
predicted.Z             0.34727    0.12546   2.768  0.00564 **
scale(RSE)             -0.18688    0.12281  -1.522  0.12807   
predicted.Z:scale(RSE)  0.11928    0.09511   1.254  0.20978   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z s(RSE)
predicted.Z -0.043              
scale(RSE)   0.012 -0.521       
prd.Z:(RSE) -0.373  0.212 -0.159
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "RSE")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(SCC) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(SCC) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4452.1   4501.1  -2218.1   4436.1     3364 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.2296 -1.0482  0.6315  0.8635  2.1392 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr 
 trait  (Intercept) 2.572e-14 1.604e-07      
 subID  (Intercept) 2.608e-01 5.107e-01      
        predicted.Z 1.120e-02 1.058e-01 -1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                        Estimate Std. Error z value Pr(>|z|)  
(Intercept)             0.297517   0.118444   2.512   0.0120 *
predicted.Z             0.279553   0.115650   2.417   0.0156 *
scale(SCC)             -0.008679   0.119455  -0.073   0.9421  
predicted.Z:scale(SCC) -0.046727   0.120280  -0.388   0.6977  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z s(SCC)
predicted.Z -0.267              
scale(SCC)   0.123 -0.318       
prd.Z:(SCC) -0.344 -0.080 -0.369
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "SCC")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(DS) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(DS) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4452.3   4501.3  -2218.2   4436.3     3364 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.2274 -1.0484  0.6330  0.8691  2.1269 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.000000 0.00000       
 subID  (Intercept) 0.259573 0.50948       
        predicted.Z 0.009853 0.09926  -1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                       Estimate Std. Error z value Pr(>|z|)  
(Intercept)            0.285057   0.116450   2.448   0.0144 *
predicted.Z            0.258527   0.115217   2.244   0.0248 *
scale(DS)             -0.020453   0.130288  -0.157   0.8753  
predicted.Z:scale(DS)  0.009681   0.117406   0.082   0.9343  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z sc(DS)
predicted.Z -0.322              
scale(DS)   -0.171  0.374       
prdc.Z:(DS)  0.313 -0.161 -0.584
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "DS")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(NFC) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(NFC) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4448.3   4497.3  -2216.1   4432.3     3364 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.2304 -1.0520  0.6456  0.8539  2.2222 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.000000 0.00000       
 subID  (Intercept) 0.233963 0.48370       
        predicted.Z 0.008746 0.09352  -1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                       Estimate Std. Error z value Pr(>|z|)   
(Intercept)             0.31018    0.10638   2.916  0.00355 **
predicted.Z             0.31111    0.11148   2.791  0.00526 **
scale(NFC)             -0.08793    0.10823  -0.812  0.41655   
predicted.Z:scale(NFC)  0.26711    0.13185   2.026  0.04278 * 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z s(NFC)
predicted.Z -0.266              
scale(NFC)  -0.041  0.035       
prd.Z:(NFC)  0.124  0.251 -0.408
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "NFC")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(SING.Ind) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(SING.Ind) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4445.3   4494.3  -2214.7   4429.3     3364 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.3660 -1.0464  0.6193  0.8716  1.9869 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr
 trait  (Intercept) 1.696e-15 4.118e-08     
 subID  (Intercept) 2.129e-01 4.614e-01     
        predicted.Z 6.052e-03 7.779e-02 1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                            Estimate Std. Error z value Pr(>|z|)   
(Intercept)                  0.22506    0.11183   2.013  0.04416 * 
predicted.Z                  0.49350    0.15106   3.267  0.00109 **
scale(SING.Ind)             -0.41608    0.15616  -2.664  0.00771 **
predicted.Z:scale(SING.Ind)  0.12930    0.06301   2.052  0.04015 * 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z s(SING
predicted.Z  0.234              
scl(SING.I) -0.124 -0.744       
p.Z:(SING.I -0.391  0.062  0.074
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "SING.Ind")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(SING.Inter) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(SING.Inter) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4450.2   4499.2  -2217.1   4434.2     3364 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.3195 -1.0424  0.6272  0.8741  2.0284 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr
 trait  (Intercept) 6.097e-15 7.808e-08     
 subID  (Intercept) 2.510e-01 5.010e-01     
        predicted.Z 1.353e-03 3.678e-02 1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                              Estimate Std. Error z value Pr(>|z|)   
(Intercept)                     0.3212     0.1107   2.900  0.00373 **
predicted.Z                     0.2784     0.1237   2.251  0.02437 * 
scale(SING.Inter)               0.1673     0.1113   1.504  0.13262   
predicted.Z:scale(SING.Inter)  -0.1145     0.0868  -1.319  0.18708   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z s(SING
predicted.Z  0.199              
scl(SING.I)  0.055 -0.040       
p.Z:(SING.I -0.205 -0.382 -0.152
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "SING.Inter")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(Proto) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(Proto) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4451.6   4500.6  -2217.8   4435.6     3364 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.2262 -1.0496  0.6326  0.8715  2.1077 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.000000 0.00000       
 subID  (Intercept) 0.246752 0.49674       
        predicted.Z 0.005939 0.07706  -1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                         Estimate Std. Error z value Pr(>|z|)  
(Intercept)               0.26086    0.11026   2.366   0.0180 *
predicted.Z               0.25533    0.11199   2.280   0.0226 *
scale(Proto)             -0.02969    0.12083  -0.246   0.8059  
predicted.Z:scale(Proto)  0.06533    0.08536   0.765   0.4441  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z scl(P)
predicted.Z -0.196              
scale(Prot)  0.093 -0.338       
prdct.Z:(P) -0.224 -0.006 -0.472
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "Proto")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(SI) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(SI) + (predicted.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4451.5   4500.4  -2217.7   4435.5     3364 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.1926 -1.0490  0.6295  0.8731  2.1288 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr 
 trait  (Intercept) 1.893e-14 1.376e-07      
 subID  (Intercept) 2.480e-01 4.980e-01      
        predicted.Z 1.495e-02 1.223e-01 -1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                      Estimate Std. Error z value Pr(>|z|)  
(Intercept)            0.27631    0.10918   2.531   0.0114 *
predicted.Z            0.22812    0.11042   2.066   0.0388 *
scale(SI)              0.08172    0.11025   0.741   0.4585  
predicted.Z:scale(SI)  0.02415    0.08217   0.294   0.7689  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z sc(SI)
predicted.Z -0.293              
scale(SI)    0.052 -0.128       
prdc.Z:(SI) -0.101 -0.344 -0.393
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "SI")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(NTB) + ( predicted.Z | subID) + ( predicted.Z | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * scale(NTB) + (predicted.Z |      subID) + (predicted.Z | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4450.7   4512.0  -2215.4   4430.7     3362 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.2315 -1.0498  0.6251  0.8688  2.0892 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr 
 trait  (Intercept) 1.100e-14 1.049e-07      
        predicted.Z 2.054e-13 4.532e-07 1.00 
 subID  (Intercept) 2.041e-01 4.518e-01      
        predicted.Z 1.935e-03 4.399e-02 -1.00
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                       Estimate Std. Error z value Pr(>|z|)   
(Intercept)             0.28922    0.09881   2.927  0.00342 **
predicted.Z             0.26267    0.10390   2.528  0.01146 * 
scale(NTB)              0.23726    0.10150   2.338  0.01941 * 
predicted.Z:scale(NTB) -0.05590    0.09787  -0.571  0.56790   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z s(NTB)
predicted.Z -0.124              
scale(NTB)   0.018 -0.049       
prd.Z:(NTB)  0.032 -0.353  0.177
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("predicted.Z", "NTB")) %>% plot()
Data were 'prettified'. Consider using `terms="predicted.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ predicted.Z*scale(RSE) + ( predicted.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(RSE) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(RSE) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4793.6   4843.2  -2388.8   4777.6     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8434 -1.0083  0.5747  0.8744  2.3956 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr 
 trait  (Intercept) 1.004e-18 1.002e-09      
 subID  (Intercept) 9.301e-01 9.644e-01      
        entropy.Z   1.594e+00 1.262e+00 -0.55
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                     Estimate Std. Error z value Pr(>|z|)   
(Intercept)            0.7127     0.2436   2.925  0.00344 **
entropy.Z             -0.5494     0.3071  -1.789  0.07362 . 
scale(RSE)            -0.2484     0.2335  -1.064  0.28751   
entropy.Z:scale(RSE)   0.3395     0.3118   1.089  0.27626   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z s(RSE)
entropy.Z   -0.412              
scale(RSE)  -0.019  0.206       
ent.Z:(RSE)  0.218 -0.047 -0.388
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "RSE")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(SCC) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(SCC) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4795.1   4844.7  -2389.6   4779.1     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8462 -1.0069  0.5800  0.8735  2.3915 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr 
 trait  (Intercept) 6.577e-20 2.564e-10      
 subID  (Intercept) 9.521e-01 9.757e-01      
        entropy.Z   1.680e+00 1.296e+00 -0.56
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                      Estimate Std. Error z value Pr(>|z|)   
(Intercept)           0.690663   0.243870   2.832  0.00462 **
entropy.Z            -0.527405   0.314479  -1.677  0.09353 . 
scale(SCC)           -0.097378   0.240594  -0.405  0.68567   
entropy.Z:scale(SCC)  0.005057   0.299650   0.017  0.98653   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z s(SCC)
entropy.Z   -0.436              
scale(SCC)  -0.071  0.122       
ent.Z:(SCC)  0.118  0.025 -0.481
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "SCC")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(DS) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(DS) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4793.8   4843.4  -2388.9   4777.8     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8406 -1.0079  0.5813  0.8733  2.4261 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.0000   0.000         
 subID  (Intercept) 0.8798   0.938         
        entropy.Z   1.6204   1.273    -0.55
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                    Estimate Std. Error z value Pr(>|z|)   
(Intercept)           0.7324     0.2379   3.078  0.00208 **
entropy.Z            -0.5417     0.3101  -1.747  0.08070 . 
scale(DS)             0.2496     0.2290   1.090  0.27583   
entropy.Z:scale(DS)  -0.3141     0.3095  -1.015  0.31015   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z sc(DS)
entropy.Z   -0.431              
scale(DS)    0.144 -0.206       
entr.Z:(DS) -0.214  0.020 -0.448
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "DS")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(NFC) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(NFC) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4791.1   4840.7  -2387.6   4775.1     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8458 -1.0084  0.5964  0.8739  2.4604 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.0000   0.0000        
 subID  (Intercept) 0.7927   0.8903        
        entropy.Z   1.2373   1.1124   -0.55
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                     Estimate Std. Error z value Pr(>|z|)   
(Intercept)            0.7115     0.2278   3.123  0.00179 **
entropy.Z             -0.5087     0.2797  -1.819  0.06897 . 
scale(NFC)             0.2597     0.2271   1.144  0.25278   
entropy.Z:scale(NFC)  -0.5659     0.2716  -2.083  0.03722 * 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z s(NFC)
entropy.Z   -0.442              
scale(NFC)   0.112 -0.175       
ent.Z:(NFC) -0.190  0.048 -0.483
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "NFC")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(SING.Ind) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(SING.Ind) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4791.7   4841.3  -2387.9   4775.7     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8319 -1.0092  0.5749  0.8759  2.4004 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.0000   0.0000        
 subID  (Intercept) 0.6746   0.8214        
        entropy.Z   1.4260   1.1941   -0.65
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                          Estimate Std. Error z value Pr(>|z|)   
(Intercept)                 0.6476     0.2177   2.974  0.00294 **
entropy.Z                  -0.6070     0.2932  -2.071  0.03840 * 
scale(SING.Ind)            -0.5182     0.2563  -2.022  0.04318 * 
entropy.Z:scale(SING.Ind)   0.1566     0.2875   0.545  0.58605   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z s(SING
entropy.Z   -0.467              
scl(SING.I) -0.005  0.244       
e.Z:(SING.I  0.228 -0.066 -0.327
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "SING.Ind")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(SING.Inter) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(SING.Inter) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4791.4   4841.0  -2387.7   4775.4     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8785 -1.0138  0.5936  0.8820  2.4011 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.0000   0.0000        
 subID  (Intercept) 0.6094   0.7806        
        entropy.Z   1.7311   1.3157   -0.51
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                            Estimate Std. Error z value Pr(>|z|)    
(Intercept)                   0.7522     0.2142   3.512 0.000444 ***
entropy.Z                    -0.4557     0.3140  -1.451 0.146682    
scale(SING.Inter)             0.5286     0.2593   2.038 0.041513 *  
entropy.Z:scale(SING.Inter)  -0.1277     0.2946  -0.433 0.664789    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z s(SING
entropy.Z   -0.382              
scl(SING.I)  0.237 -0.005       
e.Z:(SING.I -0.056  0.040 -0.301
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "SING.Inter")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(Proto) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(Proto) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4795.2   4844.8  -2389.6   4779.2     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8447 -1.0052  0.5806  0.8732  2.3944 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.0000   0.0000        
 subID  (Intercept) 0.9035   0.9505        
        entropy.Z   1.6478   1.2837   -0.54
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                       Estimate Std. Error z value Pr(>|z|)   
(Intercept)             0.69010    0.23886   2.889  0.00386 **
entropy.Z              -0.50042    0.31002  -1.614  0.10650   
scale(Proto)            0.08185    0.23908   0.342  0.73209   
entropy.Z:scale(Proto) -0.09448    0.29866  -0.316  0.75173   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z scl(P)
entropy.Z   -0.432              
scale(Prot)  0.057  0.041       
entrp.Z:(P)  0.011 -0.017 -0.281
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "Proto")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(SI) + ( entropy.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(SI) + (entropy.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4790.2   4839.8  -2387.1   4774.2     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8302 -1.0044  0.5814  0.8769  2.4447 

Random effects:
 Groups Name        Variance Std.Dev. Corr 
 trait  (Intercept) 0.0000   0.0000        
 subID  (Intercept) 0.6831   0.8265        
        entropy.Z   1.8112   1.3458   -0.64
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                    Estimate Std. Error z value Pr(>|z|)    
(Intercept)          0.74323    0.21772   3.414 0.000641 ***
entropy.Z           -0.47302    0.31535  -1.500 0.133618    
scale(SI)            0.45524    0.21902   2.079 0.037661 *  
entropy.Z:scale(SI) -0.01954    0.30810  -0.063 0.949443    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z sc(SI)
entropy.Z   -0.487              
scale(SI)    0.165 -0.075       
entr.Z:(SI) -0.024  0.071 -0.503
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "SI")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ entropy.Z*scale(NTB) + ( entropy.Z | subID) + ( entropy.Z | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ entropy.Z * scale(NTB) + (entropy.Z |      subID) + (entropy.Z | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4798.1   4860.1  -2389.0   4778.1     3630 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.8355 -1.0062  0.5801  0.8745  2.3634 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr 
 trait  (Intercept) 0.000e+00 0.000e+00      
        entropy.Z   5.571e-16 2.360e-08  NaN 
 subID  (Intercept) 8.986e-01 9.479e-01      
        entropy.Z   1.595e+00 1.263e+00 -0.57
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                     Estimate Std. Error z value Pr(>|z|)   
(Intercept)           0.68431    0.23774   2.878   0.0040 **
entropy.Z            -0.52197    0.30393  -1.717   0.0859 . 
scale(NTB)            0.24331    0.23185   1.049   0.2940   
entropy.Z:scale(NTB) -0.06185    0.29468  -0.210   0.8337   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) entr.Z s(NTB)
entropy.Z   -0.452              
scale(NTB)   0.040 -0.082       
ent.Z:(NTB) -0.066  0.031 -0.534
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("entropy.Z", "NTB")) %>% plot()
Data were 'prettified'. Consider using `terms="entropy.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(RSE) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(RSE) + (slope.Z | subID) +      (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4717.6   4767.2  -2350.8   4701.6     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.0778 -0.9845  0.4695  0.8520  2.6235 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000    0.000        
 subID  (Intercept) 8.780    2.963        
        slope.Z     8.935    2.989    0.33
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                   Estimate Std. Error z value Pr(>|z|)  
(Intercept)          0.1636     0.6400   0.256   0.7982  
slope.Z              1.2559     0.6528   1.924   0.0544 .
scale(RSE)          -1.1510     0.6472  -1.778   0.0753 .
slope.Z:scale(RSE)  -0.6067     0.6374  -0.952   0.3412  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z s(RSE)
slope.Z      0.284              
scale(RSE)  -0.004 -0.062       
slp.Z:(RSE) -0.062 -0.004  0.261
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "RSE")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(SCC) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(SCC) + (slope.Z | subID) +      (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4720.3   4769.9  -2352.2   4704.3     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.4081 -0.9850  0.4734  0.8486  2.6613 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr
 trait  (Intercept) 3.451e-16 1.858e-08     
 subID  (Intercept) 1.003e+01 3.167e+00     
        slope.Z     9.640e+00 3.105e+00 0.39
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                   Estimate Std. Error z value Pr(>|z|)  
(Intercept)          0.1138     0.6771   0.168   0.8665  
slope.Z              1.2064     0.6740   1.790   0.0735 .
scale(SCC)          -0.3103     0.6658  -0.466   0.6412  
slope.Z:scale(SCC)   0.2471     0.6620   0.373   0.7089  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z s(SCC)
slope.Z      0.342              
scale(SCC)   0.011 -0.042       
slp.Z:(SCC) -0.042  0.013  0.347
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "SCC")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(DS) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(DS) + (slope.Z | subID) +      (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4719.6   4769.2  -2351.8   4703.6     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.2079 -0.9852  0.4703  0.8481  2.6378 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000    0.000        
 subID  (Intercept) 9.636    3.104        
        slope.Z     9.323    3.053    0.34
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                  Estimate Std. Error z value Pr(>|z|)  
(Intercept)         0.1673     0.6667   0.251   0.8019  
slope.Z             1.2297     0.6656   1.848   0.0647 .
scale(DS)           0.6742     0.6705   1.006   0.3146  
slope.Z:scale(DS)   0.5111     0.6539   0.782   0.4344  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z sc(DS)
slope.Z     0.299               
scale(DS)   0.017  0.063        
slp.Z:s(DS) 0.059  0.012  0.288 
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "DS")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(NFC) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(NFC) + (slope.Z | subID) +      (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4719.1   4768.7  -2351.5   4703.1     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.3504 -0.9885  0.4707  0.8482  2.6496 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000    0.000        
 subID  (Intercept) 9.830    3.135        
        slope.Z     9.302    3.050    0.41
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                   Estimate Std. Error z value Pr(>|z|)  
(Intercept)          0.1664     0.6708   0.248   0.8041  
slope.Z              1.1872     0.6637   1.789   0.0737 .
scale(NFC)          -0.3811     0.6609  -0.577   0.5641  
slope.Z:scale(NFC)   0.6284     0.6594   0.953   0.3406  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z s(NFC)
slope.Z      0.357              
scale(NFC)  -0.017  0.010       
slp.Z:(NFC)  0.011  0.027  0.344
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "NFC")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(SING.Ind) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(SING.Ind) + (slope.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4717.4   4767.0  -2350.7   4701.4     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.1006 -0.9855  0.4721  0.8500  2.6384 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000    0.000        
 subID  (Intercept) 8.250    2.872        
        slope.Z     8.988    2.998    0.36
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                        Estimate Std. Error z value Pr(>|z|)  
(Intercept)              0.03639    0.63023   0.058   0.9540  
slope.Z                  1.33378    0.65833   2.026   0.0428 *
scale(SING.Ind)         -1.43648    0.74774  -1.921   0.0547 .
slope.Z:scale(SING.Ind) -0.27897    0.65618  -0.425   0.6707  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z s(SING
slope.Z      0.289              
scl(SING.I)  0.103 -0.132       
s.Z:(SING.I -0.124  0.028  0.091
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "SING.Ind")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(SING.Inter) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(SING.Inter) + (slope.Z |      subID) + (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4719.2   4768.8  -2351.6   4703.2     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.3813 -0.9853  0.4706  0.8482  2.6105 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000    0.000        
 subID  (Intercept) 9.447    3.074        
        slope.Z     9.766    3.125    0.37
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                          Estimate Std. Error z value Pr(>|z|)  
(Intercept)                0.23586    0.66492   0.355   0.7228  
slope.Z                    1.14832    0.67768   1.694   0.0902 .
scale(SING.Inter)          0.98467    0.75955   1.296   0.1948  
slope.Z:scale(SING.Inter)  0.08021    0.67349   0.119   0.9052  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z s(SING
slope.Z      0.317              
scl(SING.I)  0.112 -0.034       
s.Z:(SING.I -0.049  0.003  0.145
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "SING.Inter")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(Proto) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(Proto) + (slope.Z | subID) +      (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4720.7   4770.3  -2352.4   4704.7     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.3935 -0.9845  0.4706  0.8482  2.6074 

Random effects:
 Groups Name        Variance  Std.Dev.  Corr
 trait  (Intercept) 1.223e-15 3.497e-08     
 subID  (Intercept) 1.005e+01 3.170e+00     
        slope.Z     9.616e+00 3.101e+00 0.39
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                     Estimate Std. Error z value Pr(>|z|)  
(Intercept)           0.11611    0.67889   0.171   0.8642  
slope.Z               1.19302    0.67388   1.770   0.0767 .
scale(Proto)         -0.05402    0.72101  -0.075   0.9403  
slope.Z:scale(Proto)  0.23967    0.66898   0.358   0.7201  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z scl(P)
slope.Z      0.337              
scale(Prot)  0.029 -0.051       
slp.Z:sc(P) -0.053  0.000  0.270
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "Proto")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(SI) + ( slope.Z | subID) + ( 1 | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(SI) + (slope.Z | subID) +      (1 | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4718.6   4768.2  -2351.3   4702.6     3632 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.1024 -0.9805  0.4712  0.8492  2.6025 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000    0.000        
 subID  (Intercept) 9.008    3.001        
        slope.Z     9.362    3.060    0.33
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                  Estimate Std. Error z value Pr(>|z|)  
(Intercept)         0.1666     0.6465   0.258   0.7967  
slope.Z             1.1603     0.6651   1.745   0.0811 .
scale(SI)           1.0067     0.6722   1.498   0.1342  
slope.Z:scale(SI)   0.5536     0.6642   0.834   0.4045  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z sc(SI)
slope.Z      0.281              
scale(SI)    0.044 -0.025       
slp.Z:s(SI) -0.028 -0.006  0.273
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "SI")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ slope.Z*scale(NTB) + ( slope.Z | subID) + ( slope.Z | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ slope.Z * scale(NTB) + (slope.Z | subID) +      (slope.Z | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4723.3   4785.3  -2351.6   4703.3     3630 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.4383 -0.9847  0.4726  0.8491  2.5658 

Random effects:
 Groups Name        Variance Std.Dev. Corr
 trait  (Intercept) 0.000000 0.00000      
        slope.Z     0.001685 0.04104   NaN
 subID  (Intercept) 9.575622 3.09445      
        slope.Z     9.551674 3.09058  0.37
Number of obs: 3640, groups:  trait, 148; subID, 25

Fixed effects:
                   Estimate Std. Error z value Pr(>|z|)  
(Intercept)          0.1747     0.6654   0.263   0.7929  
slope.Z              1.1860     0.6712   1.767   0.0772 .
scale(NTB)           0.8431     0.6661   1.266   0.2056  
slope.Z:scale(NTB)   0.3383     0.6636   0.510   0.6102  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) slop.Z s(NTB)
slope.Z     0.327               
scale(NTB)  0.045  0.003        
slp.Z:(NTB) 0.004  0.012  0.376 
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("slope.Z", "NTB")) %>% plot()
Data were 'prettified'. Consider using `terms="slope.Z [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(RSE) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(RSE) + (scale(desirability) |  
    subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2111.4   2165.2  -1045.7   2091.4     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.8959 -1.0144  0.4983  0.9101  1.4069 

Random effects:
 Groups Name                Variance Std.Dev. Corr
 trait  (Intercept)         0.003365 0.05801      
        scale(desirability) 0.013385 0.11570  1.00
 subID  (Intercept)         0.415188 0.64435      
        scale(desirability) 0.001739 0.04170  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                               Estimate Std. Error z value Pr(>|z|)    
(Intercept)                     0.38520    0.20272   1.900  0.05741 .  
scale(desirability)             0.20932    0.06083   3.441  0.00058 ***
scale(RSE)                      0.27336    0.20314   1.346  0.17840    
scale(desirability):scale(RSE)  0.02205    0.05729   0.385  0.70036    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) s(RSE)
scl(dsrblt) 0.218               
scale(RSE)  0.019  0.019        
scl():(RSE) 0.018  0.128  0.233 
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "RSE")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(SCC) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(SCC) + (scale(desirability) |  
    subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2110.0   2163.8  -1045.0   2090.0     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.8913 -1.0078  0.4896  0.9083  1.3630 

Random effects:
 Groups Name                Variance Std.Dev. Corr
 trait  (Intercept)         0.003354 0.05791      
        scale(desirability) 0.013479 0.11610  1.00
 subID  (Intercept)         0.377675 0.61455      
        scale(desirability) 0.003107 0.05574  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                                Estimate Std. Error z value Pr(>|z|)    
(Intercept)                     0.385673   0.194130   1.987 0.046958 *  
scale(desirability)             0.209174   0.061907   3.379 0.000728 ***
scale(SCC)                      0.337654   0.195663   1.726 0.084403 .  
scale(desirability):scale(SCC) -0.006718   0.060638  -0.111 0.911783    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) s(SCC)
scl(dsrblt) 0.280               
scale(SCC)  0.026  0.024        
scl():(SCC) 0.022  0.143  0.293 
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "SCC")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(DS) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(DS) + (scale(desirability) |  
    subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2110.9   2164.7  -1045.4   2090.9     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.7639 -1.0088  0.4900  0.9098  1.4151 

Random effects:
 Groups Name                Variance  Std.Dev. Corr
 trait  (Intercept)         0.0034089 0.05839      
        scale(desirability) 0.0134907 0.11615  1.00
 subID  (Intercept)         0.3893895 0.62401      
        scale(desirability) 0.0009413 0.03068  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                              Estimate Std. Error z value Pr(>|z|)    
(Intercept)                    0.38260    0.19680   1.944 0.051891 .  
scale(desirability)            0.20787    0.06023   3.451 0.000558 ***
scale(DS)                     -0.29740    0.19610  -1.517 0.129382    
scale(desirability):scale(DS) -0.03279    0.05445  -0.602 0.547033    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) sc(DS)
scl(dsrblt)  0.167              
scale(DS)   -0.015 -0.017       
scl(d):(DS) -0.015 -0.103  0.179
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "DS")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(NFC) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(NFC) + (scale(desirability) |  
    subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2112.6   2166.4  -1046.3   2092.6     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.8577 -1.0118  0.4889  0.9094  1.3512 

Random effects:
 Groups Name                Variance Std.Dev. Corr
 trait  (Intercept)         0.003355 0.05792      
        scale(desirability) 0.013450 0.11598  1.00
 subID  (Intercept)         0.479426 0.69241      
        scale(desirability) 0.002268 0.04762  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                               Estimate Std. Error z value Pr(>|z|)    
(Intercept)                     0.38517    0.21664   1.778 0.075415 .  
scale(desirability)             0.20871    0.06126   3.407 0.000657 ***
scale(NFC)                     -0.09903    0.21536  -0.460 0.645636    
scale(desirability):scale(NFC)  0.01942    0.05410   0.359 0.719619    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) s(NFC)
scl(dsrblt)  0.244              
scale(NFC)  -0.008 -0.007       
scl():(NFC) -0.008 -0.060  0.264
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "NFC")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(SING.Ind) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(SING.Ind) +  
    (scale(desirability) | subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2112.7   2166.5  -1046.4   2092.7     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.8602 -1.0147  0.4910  0.9059  1.3673 

Random effects:
 Groups Name                Variance Std.Dev. Corr
 trait  (Intercept)         0.003337 0.05777      
        scale(desirability) 0.013431 0.11589  1.00
 subID  (Intercept)         0.486231 0.69730      
        scale(desirability) 0.002218 0.04710  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                                    Estimate Std. Error z value Pr(>|z|)    
(Intercept)                          0.38508    0.21806   1.766 0.077406 .  
scale(desirability)                  0.20921    0.06118   3.419 0.000628 ***
scale(SING.Ind)                      0.05346    0.21663   0.247 0.805083    
scale(desirability):scale(SING.Ind) -0.02489    0.05367  -0.464 0.642820    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) s(SING
scl(dsrblt) 0.242               
scl(SING.I) 0.004  0.003        
s():(SING.I 0.003  0.022  0.263 
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "SING.Ind")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(SING.Inter) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(SING.Inter) +  
    (scale(desirability) | subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2110.4   2164.2  -1045.2   2090.4     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.8256 -1.0104  0.4934  0.9125  1.4150 

Random effects:
 Groups Name                Variance  Std.Dev. Corr
 trait  (Intercept)         0.0000000 0.00000      
        scale(desirability) 0.0093185 0.09653   NaN
 subID  (Intercept)         0.3727206 0.61051      
        scale(desirability) 0.0009385 0.03063  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                                      Estimate Std. Error z value Pr(>|z|)    
(Intercept)                            0.38345    0.19280   1.989 0.046716 *  
scale(desirability)                    0.20642    0.05956   3.466 0.000529 ***
scale(SING.Inter)                     -0.33236    0.19353  -1.717 0.085910 .  
scale(desirability):scale(SING.Inter) -0.03516    0.05701  -0.617 0.537466    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) s(SING
scl(dsrblt)  0.162              
scl(SING.I) -0.022 -0.021       
s():(SING.I -0.020 -0.138  0.177
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "SING.Inter")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(Proto) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(Proto) +  
    (scale(desirability) | subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2109.2   2163.0  -1044.6   2089.2     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.3834 -1.0152  0.5623  0.9062  1.4029 

Random effects:
 Groups Name                Variance Std.Dev. Corr
 trait  (Intercept)         0.003279 0.05726      
        scale(desirability) 0.013172 0.11477  1.00
 subID  (Intercept)         0.374479 0.61195      
        scale(desirability) 0.001273 0.03567  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                                 Estimate Std. Error z value Pr(>|z|)    
(Intercept)                       0.39167    0.19357   2.023 0.043030 *  
scale(desirability)               0.21533    0.06089   3.536 0.000406 ***
scale(Proto)                     -0.38362    0.20188  -1.900 0.057402 .  
scale(desirability):scale(Proto) -0.07237    0.07165  -1.010 0.312464    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) scl(P)
scl(dsrblt)  0.194              
scale(Prot) -0.043 -0.047       
scl(ds):(P) -0.038 -0.225  0.229
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "Proto")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(SI) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(SI) + (scale(desirability) |  
    subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2110.7   2164.5  -1045.4   2090.7     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-4.0641 -1.0169  0.5313  0.9074  1.5565 

Random effects:
 Groups Name                Variance Std.Dev. Corr
 trait  (Intercept)         0.003313 0.05756      
        scale(desirability) 0.013506 0.11621  1.00
 subID  (Intercept)         0.492358 0.70168      
        scale(desirability) 0.002593 0.05092  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                              Estimate Std. Error z value Pr(>|z|)    
(Intercept)                    0.38629    0.21938   1.761 0.078272 .  
scale(desirability)            0.21359    0.06156   3.470 0.000521 ***
scale(SI)                     -0.05538    0.21816  -0.254 0.799632    
scale(desirability):scale(SI) -0.08211    0.05464  -1.503 0.132865    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) sc(SI)
scl(dsrblt)  0.260              
scale(SI)   -0.008 -0.009       
scl(d):(SI) -0.009 -0.079  0.280
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "SI")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

m <- glmer( as.factor(ingChoiceN) ~ scale(desirability)*scale(NTB) + ( scale(desirability) | subID) + ( scale(desirability) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
boundary (singular) fit: see help('isSingular')
summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ scale(desirability) * scale(NTB) + (scale(desirability) |  
    subID) + (scale(desirability) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  2112.8   2166.6  -1046.4   2092.8     1594 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-3.9039 -1.0131  0.5005  0.9074  1.4211 

Random effects:
 Groups Name                Variance Std.Dev. Corr
 trait  (Intercept)         0.00337  0.05805      
        scale(desirability) 0.01338  0.11569  1.00
 subID  (Intercept)         0.47789  0.69130      
        scale(desirability) 0.00211  0.04593  1.00
Number of obs: 1604, groups:  trait, 148; subID, 11

Fixed effects:
                               Estimate Std. Error z value Pr(>|z|)    
(Intercept)                     0.38539    0.21633   1.782  0.07483 .  
scale(desirability)             0.20947    0.06112   3.427  0.00061 ***
scale(NTB)                     -0.11353    0.21542  -0.527  0.59819    
scale(desirability):scale(NTB) -0.01384    0.05480  -0.253  0.80060    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) scl(d) s(NTB)
scl(dsrblt)  0.237              
scale(NTB)  -0.010 -0.010       
scl():(NTB) -0.009 -0.075  0.255
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
ggpredict(m, c("desirability", "NTB")) %>% plot()
Data were 'prettified'. Consider using `terms="desirability [all]"` to get smooth plots.

summary(m)
Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
 Family: binomial  ( logit )
Formula: as.factor(ingChoiceN) ~ predicted.Z * novel * scale(NFC) + (predicted.Z +  
    novel | subID) + (SE.Z * as.factor(novel) | trait)
   Data: fullTest
Control: glmerControl(optimizer = "bobyqa", optCtrl = list(maxfun = 1e+05))

     AIC      BIC   logLik deviance df.resid 
  4476.7   4623.6  -2214.3   4428.7     3348 

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.2596 -1.0499  0.6322  0.8416  2.2666 

Random effects:
 Groups Name                          Variance Std.Dev. Corr             
 trait  (Intercept)                   0.000000 0.00000                   
        SE.Z                          0.001061 0.03257    NaN            
        as.factor(novel)Held Out      0.002307 0.04803    NaN -1.00      
        SE.Z:as.factor(novel)Held Out 0.038538 0.19631    NaN -1.00  1.00
 subID  (Intercept)                   0.257661 0.50760                   
        predicted.Z                   0.008917 0.09443  -1.00            
        novelHeld Out                 0.002718 0.05214  -1.00  1.00      
Number of obs: 3372, groups:  trait, 148; subID, 25

Fixed effects:
                                      Estimate Std. Error z value Pr(>|z|)   
(Intercept)                           0.298628   0.116464   2.564  0.01034 * 
predicted.Z                           0.318133   0.123109   2.584  0.00976 **
novelHeld Out                         0.031633   0.078956   0.401  0.68869   
scale(NFC)                           -0.101568   0.117287  -0.866  0.38650   
predicted.Z:novelHeld Out            -0.007964   0.083210  -0.096  0.92376   
predicted.Z:scale(NFC)                0.225593   0.143761   1.569  0.11660   
novelHeld Out:scale(NFC)              0.037730   0.075930   0.497  0.61925   
predicted.Z:novelHeld Out:scale(NFC)  0.108722   0.096644   1.125  0.26060   
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) prdc.Z nvlHlO s(NFC) pr.Z:HO p.Z:(N nHO:(N
predicted.Z -0.260                                           
novelHeldOt -0.409  0.048                                    
scale(NFC)  -0.018  0.017 -0.036                             
prdctd.Z:HO  0.057 -0.424 -0.001  0.025                      
prd.Z:(NFC)  0.097  0.270 -0.001 -0.400 -0.138               
nvlHO:(NFC) -0.040  0.027  0.084 -0.393 -0.037   0.129       
p.Z:HO:(NFC  0.003 -0.139  0.046  0.124  0.312  -0.405 -0.259
optimizer (bobyqa) convergence code: 0 (OK)
boundary (singular) fit: see help('isSingular')
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(SCC) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE","novel","SCC")) %>% plot()
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(DS) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE","novel","DS")) %>% plot()
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(NFC) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE","novel", "NFC")) %>% plot()
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(SING.Ind) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE","novel","SING.Ind")) %>% plot()
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(SING.Inter) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE","novel","SING.Inter")) %>% plot()
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(Proto) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE", "Proto")) %>% plot()
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(SI) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE","novel","SI")) %>% plot()
m <- glmer( as.factor(ingChoiceN) ~ SE.Z*as.factor(novel)*scale(NTB) + ( SE.Z+as.factor(novel) | subID) + ( SE.Z+as.factor(novel) | trait), data = fullTest, family = binomial, control = glmerControl(optimizer = "bobyqa",
                                    optCtrl = list(maxfun = 100000)),
    nAGQ = 1)
summary(m)
ggpredict(m, c("SE", "novel", "NTB")) %>% plot()
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkobG1lclRlc3QpCmxpYnJhcnkoZ2dlZmZlY3RzKQpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KHJlcG9ydCkKbGlicmFyeShyMmdsbW0pCmxpYnJhcnkocGVyZm9ybWFuY2UpCmxpYnJhcnkoYnJvb20ubWl4ZWQpCmBgYAoKYGBge3J9CmZ1bGxUZXN0IDwtIHJlYWQuY3N2KCIuLi9DbGVhbmluZy9vdXRwdXQvZnVsbFRlc3QuY3N2IikKZnVsbFRyYWluIDwtIHJlYWQuY3N2KCIuLi9DbGVhbmluZy9vdXRwdXQvZnVsbFRyYWluLmNzdiIpCgp0cmFpdHNGcmVxcyA8LSByZWFkLmNzdigiLi4vQ2xlYW5pbmcvb3V0cHV0L3RyYWl0RnJlcU92ZXJVbmRlci5jc3YiKQp0cmFpdHNGcmVxcyA8LSByZW5hbWUodHJhaXRzRnJlcXMsIHByb3BzID0gb3B0aW9uQ2hvaWNlTikKZnVsbFRlc3QgPC0gbWVyZ2UoZnVsbFRlc3QsIHRyYWl0c0ZyZXFzW2MoInRyYWl0IiwicHJvcHMiKV0sIGJ5ID0gInRyYWl0IikKZnVsbFRlc3QkcHJvcENvcnIgPC0gaWZlbHNlKGZ1bGxUZXN0JEVzdGltYXRvcj09IlVuZGVyZXN0aW1hdG9yIiwgMS1mdWxsVGVzdCRwcm9wcywgZnVsbFRlc3QkcHJvcHMpCgp1U3VicyA8LSB1bmlxdWUoZnVsbFRlc3Qkc3ViSUQpCgppbmREaWZmcyA8LSBmdWxsVGVzdFshZHVwbGljYXRlZChmdWxsVGVzdCRzdWJJRCksXQpgYGAKCmBgYHtyfQpmdWxsVGVzdCRpbmdDaG9pY2VOIDwtIGFzLmZhY3RvcihmdWxsVGVzdCRpbmdDaG9pY2VOKQpmdWxsVGVzdCRub3ZlbCA8LSBhcy5mYWN0b3IoZnVsbFRlc3Qkbm92ZWwpCmZ1bGxUZXN0JHNlbGZSZXNwLlogPC0gc2NhbGUoZnVsbFRlc3Qkc2VsZlJlc3ApCmZ1bGxUZXN0JFNFLlogPC0gc2NhbGUoZnVsbFRlc3QkU0UpCmZ1bGxUZXN0JGlTRS5aIDwtIHNjYWxlKGZ1bGxUZXN0JGlTRSkKZnVsbFRlc3Qkb1NFLlogPC0gc2NhbGUoZnVsbFRlc3Qkb1NFKQpmdWxsVGVzdCRwcmVkaWN0ZWQuWiA8LSBzY2FsZShmdWxsVGVzdCRwcmVkaWN0ZWQpCmZ1bGxUZXN0JHNsb3BlLlogPC0gc2NhbGUoZnVsbFRlc3Qkc2xvcGUpCmZ1bGxUZXN0JGVudHJvcHkuWiA8LSBzY2FsZShmdWxsVGVzdCRlbnRyb3B5KQpmdWxsVGVzdCRXU1IuWiA8LSBzY2FsZShmdWxsVGVzdCRXU1IpCmZ1bGxUZXN0JG5laWdoQXZlT3V0U0UuWiA8LSBzY2FsZShmdWxsVGVzdCRuZWlnaEF2ZU91dFNFKQpmdWxsVGVzdCRuZWlnaEF2ZUFsbFNFLlogPC0gc2NhbGUoZnVsbFRlc3QkbmVpZ2hBdmVBbGxTRSkKZnVsbFRlc3QkbmVpZ2hBdmVJblNFLlogPC0gc2NhbGUoZnVsbFRlc3QkbmVpZ2hBdmVJblNFKQpgYGAKCgpgYGB7cn0KZnVsbFRlc3Qkbm92ZWwgPC0gYXMuZmFjdG9yKGZ1bGxUZXN0JG5vdmVsKQpsZXZlbHMoZnVsbFRlc3Qkbm92ZWwpIDwtIGxpc3QoIlRyYWluZWQiICA9ICIwIiwgIkhlbGQgT3V0IiA9ICIxIikKYGBgCgoKYGBge3J9CiMgcHJvcC50ZXN0KHRyYWl0c0ZyZXFzJG9wdGlvbkNob2ljZU4sIHRyYWl0c0ZyZXFzJE4sIHA9cmVwKC41LGxlbmd0aCh0cmFpdHNGcmVxcyROKSkpCiMgCiMgbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IHRyYWl0ICsgKCAxIHwgc3ViSUQpICwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKIyAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKIyAgICAgbkFHUSA9IDEpCiMgCiMgZnVsbFRlc3QkdHJhaXQgPC0gYXMuZmFjdG9yKGZ1bGxUZXN0JHRyYWl0KQojIGNvbnRyYXN0cyhmdWxsVGVzdCR0cmFpdCkgPC0gY29udHIuc3VtKDE0OCkKIyBtIDwtIGdsbShvcHRpb25DaG9pY2VOIH4gdHJhaXQsIGZhbWlseSA9IGJpbm9taWFsLAojICAgICAgICAgICBkYXRhID0gZnVsbFRlc3QKIyAgICAgICAgICApCiMgc3VtbWFyeShtKQpgYGAKCgpgYGB7cn0KIyB0cmFpdHNGcmVxcyR0cmFpdCA8LSBhcy5mYWN0b3IodHJhaXRzRnJlcXMkdHJhaXQpCiMgY29udHJhc3RzKHRyYWl0c0ZyZXFzJHRyYWl0KSA8LSBjb250ci5zdW0oMTQ4KQojIG0gPC0gZ2xtKG9wdGlvbkNob2ljZU4gfiAxLCBmYW1pbHkgPSBiaW5vbWlhbCwKIyAgICAgICAgICAgZGF0YSA9IHRyYWl0c0ZyZXFzCiMgICAgICAgICAgKQojIHN1bW1hcnkobSkKIyBjaGVja19vdmVyZGlzcGVyc2lvbihtKQojIGNoZWNrX21vZGVsKG0pCiMgCiMgbSA8LSBnbG0ob3B0aW9uQ2hvaWNlTiB+IHRyYWl0LCBmYW1pbHkgPSBxdWFzaWJpbm9taWFsLAojICAgICAgICAgICBkYXRhID0gdHJhaXRzRnJlcXMKIyAgICAgICAgICApCiMgY2hlY2tfb3ZlcmRpc3BlcnNpb24obSkKIyBjaGVja19tb2RlbChtKQojIAojIHQudGVzdCgpCiMgCiMgbSA8LSBnbG0ob3B0aW9uQ2hvaWNlTiB+IHRyYWl0LCBmYW1pbHkgPSBwb2lzc29uLAojICAgICAgICAgICBkYXRhID0gdHJhaXRzRnJlcXMKIyAgICAgICAgICApCiMgY2hlY2tfb3ZlcmRpc3BlcnNpb24obSkKIyBjaGVja19tb2RlbChtKQpgYGAKCmBgYHtyfQpwcm9wTWF0cml4IDwtIG1hdHJpeChucm93PTE0OCxuY29sPTcpCmZvcihpIGluIDE6MTQ4KXsKICAgIHRyYWl0RGYgPC0gc3Vic2V0KGZ1bGxUZXN0LCBJZHg9PWkpCiAgICB0ZXN0IDwtIHQudGVzdChhcy5udW1lcmljKHRyYWl0RGYkaW5nQ2hvaWNlTiktMSwgbXU9LjUwKQogICAgcHJvcE1hdHJpeFtpLCBdIDwtIGMoaSwgdGVzdCRzdGF0aXN0aWMsIHRlc3QkcC52YWx1ZSwgdGVzdCRjb25mLmludCwgdGVzdCRlc3RpbWF0ZSwgdGVzdCRwYXJhbWV0ZXIpCn0KY29sbmFtZXMocHJvcE1hdHJpeCkgPC0gYygiSWR4IiwgInN0YXQiLCAicCIsICJMQ0kiLCAiVUNJIiwgImVzdCIsICJwYXJhbSIpCnByb3BNYXRyaXggPC0gYXMuZGF0YS5mcmFtZShwcm9wTWF0cml4KQpwcm9wTWF0cml4JHRyYWl0IDwtIHRyYWl0c0ZyZXFzJHRyYWl0CnByb3BNYXRyaXggPC0gcHJvcE1hdHJpeFtvcmRlcihwcm9wTWF0cml4JHApLF0KYGBgCgpgYGB7cn0KbGlicmFyeShjb3JycikKeCA8LSBpbmREaWZmcyAlPiUgCiAgICBzZWxlY3QoZ3JvdXBIb21vcGgsIERTOlNJTkcuSW50ZXIpICU+JQogIGNvcnJlbGF0ZSgpICU+JSAKICBmb2N1cyhncm91cEhvbW9waCkgJT4lCiAgICBhcnJhbmdlKGdyb3VwSG9tb3BoKQoKeCAlPiUgCiAgbXV0YXRlKHJvd25hbWUgPSBmYWN0b3Iocm93bmFtZSwgbGV2ZWxzID0gcm93bmFtZVtvcmRlcihncm91cEhvbW9waCldKSkgJT4lICAjIE9yZGVyIGJ5IGNvcnJlbGF0aW9uIHN0cmVuZ3RoCiAgZ2dwbG90KGFlcyh4ID0gcm93bmFtZSwgeSA9IGdyb3VwSG9tb3BoKSkgKwogICAgZ2VvbV9iYXIoc3RhdCA9ICJpZGVudGl0eSIpICsKICAgIHlsYWIoIkNvcnJlbGF0aW9uIENvZWZmaWNpZW50IikgKwogICAgeGxhYigiSW5kaXZpZHVhbCBEaWZmZXJlbmNlcyIpICsgdGhlbWVfZ3JleShiYXNlX3NpemUgPSA5KSAgKyB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChhbmdsZSA9IDkwLGhqdXN0ID0gMSkpICsKICB0aGVtZShheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dCggCiAgICAgICAgICAgICAgICAgICAgICAgICAgIHNpemUgPSA5LCBhbmdsZSA9IDQ1LCB2anVzdCA9IDEpKSArIHRoZW1lKGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dCh2anVzdD0xLjkpKSArIHRoZW1lKGF4aXMudGV4dD1lbGVtZW50X3RleHQoc2l6ZT05KSwKICAgICAgICBheGlzLnRpdGxlPWVsZW1lbnRfdGV4dChzaXplPTksZmFjZT0iYm9sZCIpKSArIHRoZW1lKGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9OSkpICsgdGhlbWUocGFuZWwuYm9yZGVyID0gZWxlbWVudF9yZWN0KGNvbG91ciA9ICJibGFjayIsIGZpbGwgPSBOQSwgc2l6ZSA9MSkpICsgdGhlbWUobGVnZW5kLnRpdGxlID0gZWxlbWVudF9ibGFuaygpKSArIHRoZW1lKHBhbmVsLmdyaWQubWFqb3IgPSBlbGVtZW50X2JsYW5rKCksIHBhbmVsLmdyaWQubWlub3IgPSBlbGVtZW50X2JsYW5rKCksCnBhbmVsLmJhY2tncm91bmQgPSBlbGVtZW50X2JsYW5rKCksIGF4aXMubGluZSA9IGVsZW1lbnRfbGluZShjb2xvdXIgPSAiYmxhY2siKSkKYGBgCiMgRmFtaWxpYXJpdHkgcHJlZGljdHMgUmVhY3Rpb24gVGltZQoKYGBge3J9Cm0gPC0gbG1lciggbG9nKFJUKSB+IGZhbSArICggZmFtIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoImZhbSIpKSAlPiUgcGxvdCgpCnRpZHkobSxjb25mLmludD1UUlVFLGV4cG9uZW50aWF0ZT1UUlVFLGVmZmVjdHM9ImZpeGVkIikKYGBgCgpgYGB7cn0KbSA8LSBsbWVyKCBsb2coUlQpIH4gZmFtICsgcHJvcENvcnIgKyBkZXNpcmFiaWxpdHkgKyAoIGZhbSB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0KQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJmYW0iKSkgJT4lIHBsb3QoKQp0aWR5KG0sY29uZi5pbnQ9VFJVRSxleHBvbmVudGlhdGU9VFJVRSxlZmZlY3RzPSJmaXhlZCIpCmBgYAoKCiMgUmVwbGljYXRpb24gb2YgcHJpb3Igc2VsZi1hbmNob3JpbmcgZmluZGluZ3M6IFNlbGYtZXZhbHVhdGlvbnMgcHJlZGljdGluZyBpbmdyb3VwIGV2YWx1YXRpb25zCgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IHNlbGZSZXNwLlogKyBwcm9wQ29yciArIGRlc2lyYWJpbGl0eSArICggc2VsZlJlc3AuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpsaWJyYXJ5KGJyb29tLm1peGVkKQp0aWR5KG0sY29uZi5pbnQ9VFJVRSxleHBvbmVudGlhdGU9VFJVRSxlZmZlY3RzPSJmaXhlZCIpCnIyYmV0YShtKQpnZ3ByZWRpY3QobSwgYygic2VsZlJlc3AuWiIpKSAlPiUgcGxvdChzaG93LnRpdGxlPUYpICsgeGxhYigiU2VsZi1EZXNjcmlwdGl2ZW5lc3MiKSArIHlsYWIoIkxpa2VsaWhvb2Qgb2YgSW5ncm91cCBDaG9pY2UiKSArIGp0b29sczo6dGhlbWVfYXBhKCkKZ2dzYXZlKCJ+L0RvY3VtZW50cy9VQyBSaXZlcnNpZGUvU3R1ZGllcy9TZWxmLUFuY2hvcmluZy9GaWd1cmVzL1NlbGZQcm9qZWN0aW9uLnRpZmYiLGRwaT02MDApCmBgYAoKIyBEb2VzIHNpbWlsYXJpdHktd2VpZ2h0ZWQgc2VsZi1ldmFsdWF0aW9uIGF2ZXJhZ2UgcHJlZGljdCBpbmdyb3VwIGNob2ljZXM/CgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IFdTUi5aICsgcHJvcENvcnIgKyBkZXNpcmFiaWxpdHkgKyAoIFdTUi5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCnRpZHkobSxjb25mLmludD1UUlVFLGV4cG9uZW50aWF0ZT1UUlVFLGVmZmVjdHM9ImZpeGVkIikKcjJiZXRhKG0pCmBgYAojIERvZXMgc2VsZi1ldmFsdWF0aW9uIHdlaWdodGVkIHNpbWlsYXJpdHkgcHJlZGljdCBpbmdyb3VwIGNob2ljZXM/CgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IFNFLlogKyBwcm9wQ29yciArIGRlc2lyYWJpbGl0eSArICggU0UuWiB8IHN1YklEKSAsIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKdGlkeShtLGNvbmYuaW50PVRSVUUsZXhwb25lbnRpYXRlPVRSVUUsZWZmZWN0cz0iZml4ZWQiKQpyMmJldGEobSkKYGBgCgojIERvIGNyb3NzLXZhbGlkYXRlZCBzaW1pbGFyaXR5KnNlbGYtZXZhbHVhdGlvbiBwcmVkaWN0aW9ucyBwcmVkaWN0IGluZ3JvdXAgY2hvaWNlcz8gCgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IHByZWRpY3RlZC5aICsgcHJvcENvcnIgKyBkZXNpcmFiaWxpdHkgKyAoIHByZWRpY3RlZC5aIHwgc3ViSUQpICsgKDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKdGlkeShtLGNvbmYuaW50PVRSVUUsZXhwb25lbnRpYXRlPVRSVUUsZWZmZWN0cz0iZml4ZWQiKQpyMmJldGEobSkKZ2dwcmVkaWN0KG0sIGMoInByZWRpY3RlZC5aIikpICU+JSBwbG90KHNob3cudGl0bGU9RikgKyB4bGFiKCJDcm9zcy1WYWxpZGF0ZWQgU2VsZi1EZXNjcmlwdGl2ZW5lc3MgUHJlZGljdGlvbnMiKSArIHlsYWIoIkxpa2VsaWhvb2Qgb2YgSW5ncm91cCBDaG9pY2UiKSArIGp0b29sczo6dGhlbWVfYXBhKCkKZ2dzYXZlKCJ+L0RvY3VtZW50cy9VQyBSaXZlcnNpZGUvU3R1ZGllcy9TZWxmLUFuY2hvcmluZy9GaWd1cmVzL0NWcHJlZGljdGlvbi50aWZmIixkcGk9NjAwKQpgYGAKCiMgRG8gY3Jvc3MtdmFsaWRhdGVkIHNpbWlsYXJpdHkqc2VsZi1ldmFsdWF0aW9uIHByZWRpY3Rpb25zIHByZWRpY3QgaW5ncm91cCBjaG9pY2VzLCByZWdhcmRsZXNzIG9mIHdoZXRoZXIgaXQgd2FzIHNlZW4gcHJpb3Igb3Igbm90PwoKYGBge3J9Cm0gPC0gZ2xtZXIoIGluZ0Nob2ljZU4gfiBwcmVkaWN0ZWQuWiAqIG5vdmVsICsgKCBwcmVkaWN0ZWQuWiArIG5vdmVsIHwgc3ViSUQpICwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQp0aWR5KG0sY29uZi5pbnQ9VFJVRSxleHBvbmVudGlhdGU9VFJVRSxlZmZlY3RzPSJmaXhlZCIpCnIyYmV0YShtKQpnZ3ByZWRpY3QobSwgYygicHJlZGljdGVkLloiLCAibm92ZWwiKSkgJT4lIHBsb3Qoc2hvdy50aXRsZT1GKSsgeGxhYigiQ3Jvc3MtVmFsaWRhdGVkIFNlbGYtRGVzY3JpcHRpdmVuZXNzIikgKyB5bGFiKCJMaWtlbGlob29kIG9mIEluZ3JvdXAgQ2hvaWNlIikgKyBqdG9vbHM6OnRoZW1lX2FwYSgpICArIHNjYWxlX2NvbG9yX2Rpc2NyZXRlKGxhYmVscyA9IGMoIlRyYWluZWQiLCJIZWxkLU91dCIpKQpnZ3NhdmUoIn4vRG9jdW1lbnRzL1VDIFJpdmVyc2lkZS9TdHVkaWVzL1NlbGYtQW5jaG9yaW5nL0ZpZ3VyZXMvQ1ZwcmVkaWN0aW9uR2VuZXJhbGl6ZS50aWZmIixkcGk9NjAwKQpgYGAKCiMgRG9lcyBnZW5lcmFsaXphdGlvbiBkZXBlbmQgb24gb3V0ZGVncmVlPwoKYGBge3J9Cm0gPC0gZ2xtZXIoIGluZ0Nob2ljZU4gfiBwcmVkaWN0ZWQuWiAqIG5vdmVsICogb3V0RGVncmVlICArICggcHJlZGljdGVkLlogKyBub3ZlbCB8IHN1YklEKSAsIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKdGlkeShtLGNvbmYuaW50PVRSVUUsZXhwb25lbnRpYXRlPVRSVUUsZWZmZWN0cz0iZml4ZWQiKQpyMmJldGEobSkKZ2dwcmVkaWN0KG0sIGMoInByZWRpY3RlZC5aIiwgIm91dERlZ3JlZSIgLCJub3ZlbCIpKSAlPiUgcGxvdChzaG93LnRpdGxlPUYpKyB4bGFiKCJDcm9zcy1WYWxpZGF0ZWQgU2VsZi1EZXNjcmlwdGl2ZW5lc3MiKSArIHlsYWIoIkxpa2VsaWhvb2Qgb2YgSW5ncm91cCBDaG9pY2UiKSArIGp0b29sczo6dGhlbWVfYXBhKCkKZ2dzYXZlKCJ+L0RvY3VtZW50cy9VQyBSaXZlcnNpZGUvU3R1ZGllcy9TZWxmLUFuY2hvcmluZy9GaWd1cmVzL0NWcHJlZGljdGlvbkdlbmVyYWxpemUudGlmZiIsZHBpPTYwMCkKYGBgCiMgRG9lcyBnZW5lcmFsaXphdGlvbiBkZXBlbmQgb24gaW5kZWdyZWU/CgpOb3BlCgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IHByZWRpY3RlZC5aICogbm92ZWwgKiBpbkRlZ3JlZSArICggcHJlZGljdGVkLlogKyBub3ZlbCB8IHN1YklEKSAsIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKdGlkeShtLGNvbmYuaW50PVRSVUUsZXhwb25lbnRpYXRlPVRSVUUsZWZmZWN0cz0iZml4ZWQiKQpyMmJldGEobSkKZ2dwcmVkaWN0KG0sIGMoInByZWRpY3RlZC5aIiwgIm91dERlZ3JlZSIgLCJub3ZlbCIpKSAlPiUgcGxvdChzaG93LnRpdGxlPUYpKyB4bGFiKCJDcm9zcy1WYWxpZGF0ZWQgU2VsZi1EZXNjcmlwdGl2ZW5lc3MiKSArIHlsYWIoIkxpa2VsaWhvb2Qgb2YgSW5ncm91cCBDaG9pY2UiKSArIGp0b29sczo6dGhlbWVfYXBhKCkKI2dnc2F2ZSgifi9Eb2N1bWVudHMvVUMgUml2ZXJzaWRlL1N0dWRpZXMvU2VsZi1BbmNob3JpbmcvRmlndXJlcy9DVnByZWRpY3Rpb25HZW5lcmFsaXplLnRpZmYiLGRwaT02MDApCmBgYAoKIyBOZWlnaGJvcmluZyBEZXBlbmRlbmNpZXMgUHJlZGljdGluZyBDaG9pY2VzCgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IG5laWdoQXZlT3V0U0UuWiAgKyBwcm9wQ29yciArIGRlc2lyYWJpbGl0eSArICggbmVpZ2hBdmVPdXRTRS5aIHwgc3ViSUQpICsgKDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKdGlkeShtLGNvbmYuaW50PVRSVUUsZXhwb25lbnRpYXRlPVRSVUUsZWZmZWN0cz0iZml4ZWQiKQpyMmJldGEobSkKZ2dwcmVkaWN0KG0sIGMoIm5laWdoQXZlT3V0U0UuWiIpKSAlPiUgcGxvdChzaG93LnRpdGxlPUYpKyB4bGFiKCJPdXR3YXJkcyBOZWlnaGJvcmluZyBTZWxmLUV2YWx1YXRpb25zIikgKyB5bGFiKCJMaWtlbGlob29kIG9mIEluZ3JvdXAgQ2hvaWNlIikgKyBqdG9vbHM6OnRoZW1lX2FwYSgpCmdnc2F2ZSgifi9Eb2N1bWVudHMvVUMgUml2ZXJzaWRlL1N0dWRpZXMvU2VsZi1BbmNob3JpbmcvRmlndXJlcy9vdXRkZWdyZWVOZWlnaGJvcnMudGlmZiIsZHBpPTYwMCkKYGBgCgojIEdlbmVyYWxpemF0aW9uIG9mIE91dGRlZ3JlZSBOZWlnaGJvcmluZyBTZWxmLUV2YWx1YXRpb25zCgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IG5laWdoQXZlT3V0U0UuWiAqIG5vdmVsICsgKCBuZWlnaEF2ZU91dFNFLlogKyBub3ZlbCB8IHN1YklEKSArICgxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCnRpZHkobSxjb25mLmludD1UUlVFLGV4cG9uZW50aWF0ZT1UUlVFLGVmZmVjdHM9ImZpeGVkIikKcjJiZXRhKG0pCmdncHJlZGljdChtLCBjKCJuZWlnaEF2ZU91dFNFLloiLCJub3ZlbCIpKSAlPiUgcGxvdChzaG93LnRpdGxlPUYpKyB4bGFiKCJPdXR3YXJkcyBOZWlnaGJvcmluZyBTZWxmLUV2YWx1YXRpb25zIikgKyB5bGFiKCJMaWtlbGlob29kIG9mIEluZ3JvdXAgQ2hvaWNlIikgKyBqdG9vbHM6OnRoZW1lX2FwYSgpCmdnc2F2ZSgifi9Eb2N1bWVudHMvVUMgUml2ZXJzaWRlL1N0dWRpZXMvU2VsZi1BbmNob3JpbmcvRmlndXJlcy9vdXRkZWdyZWVOZWlnaGJvcnNHZW5lcmFsaXphdGlvbi50aWZmIixkcGk9NjAwKQpgYGAKCiMgRG9lcyBlbnRyb3B5IChpLmUuLCB1bmNlcnRhaW50eSkgcHJlZGljdCBsaWtlbGlob29kIG9mIGluZ3JvdXAgY2hvaWNlcz8KCmBgYHtyfQptIDwtIGdsbWVyKCBpbmdDaG9pY2VOIH4gZW50cm9weS5aICArIHByb3BDb3JyICsgZGVzaXJhYmlsaXR5ICsgKCBlbnRyb3B5LlogfCBzdWJJRCkgKyAoMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQp0aWR5KG0sY29uZi5pbnQ9VFJVRSxleHBvbmVudGlhdGU9VFJVRSxlZmZlY3RzPSJmaXhlZCIpCnIyYmV0YShtKQpnZ3ByZWRpY3QobSwgYygiZW50cm9weS5aIikpICU+JSBwbG90KHNob3cudGl0bGU9RikgKyB4bGFiKCJVbmNlcnRhaW50eSIpICsgeWxhYigiTGlrZWxpaG9vZCBvZiBJbmdyb3VwIENob2ljZSIpICsganRvb2xzOjp0aGVtZV9hcGEoKQpnZ3NhdmUoIn4vRG9jdW1lbnRzL1VDIFJpdmVyc2lkZS9TdHVkaWVzL1NlbGYtQW5jaG9yaW5nL0ZpZ3VyZXMvVW5jZXJ0YWludHkudGlmZiIsZHBpPTYwMCkKYGBgCgojIERvZXMgYSBsaW5lYXIgdHJlbmQgb2Ygc2ltaWxhcml0eS1iYXNlZCBwcm9iYWJpbGl0aWVzIHByZWRpY3QgaW5ncm91cCBjaG9pY2VzPwoKYGBge3J9Cm0gPC0gZ2xtZXIoIGluZ0Nob2ljZU4gfiBzbG9wZS5aICArIHByb3BDb3JyICsgZGVzaXJhYmlsaXR5ICsgKCBzbG9wZS5aIHwgc3ViSUQpICsgKDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKdGlkeShtLGNvbmYuaW50PVRSVUUsZXhwb25lbnRpYXRlPVRSVUUsZWZmZWN0cz0iZml4ZWQiKQpnZ3ByZWRpY3QobSwgYygic2xvcGUuWiIpKSAlPiUgcGxvdChzaG93LnRpdGxlPUYpICsgeGxhYigiTGluZWFyIFRyZW5kIG9mIEdyZWF0ZXIgU2VsZi1EZXNjcmlwdGl2ZW5lc3MiKSArIHlsYWIoIkxpa2VsaWhvb2Qgb2YgSW5ncm91cCBDaG9pY2UiKSArIGp0b29sczo6dGhlbWVfYXBhKCkKZ2dzYXZlKCJ+L0RvY3VtZW50cy9VQyBSaXZlcnNpZGUvU3R1ZGllcy9TZWxmLUFuY2hvcmluZy9GaWd1cmVzL1Nsb3BlLnRpZmYiLGRwaT02MDApCmBgYAoKIyBEb2VzIGEgbGluZWFyIHRyZW5kIG9mIHNpbWlsYXJpdHktYmFzZWQgcHJvYmFiaWxpdGllcyBwcmVkaWN0IGluZ3JvdXAgY2hvaWNlcywgY29udHJvbGxpbmcgZm9yIHNlbGYtZGVzY3JpcHRpdmVuZXNzPwoKYGBge3J9Cm0gPC0gZ2xtZXIoIGluZ0Nob2ljZU4gfiBzY2FsZShzbG9wZSkgKyBzZWxmUmVzcC5aICArICggc2NhbGUoc2xvcGUpICsgc2VsZlJlc3AuWiB8IHN1YklEKSArICgxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCnRpZHkobSxjb25mLmludD1UUlVFLGV4cG9uZW50aWF0ZT1UUlVFLGVmZmVjdHM9ImZpeGVkIikKYGBgCgojIERvZXMgYSBub24tcGFyYW1ldHJpYyB0cmVuZCBvZiBzaW1pbGFyaXR5LWJhc2VkIHByb2JhYmlsaXRpZXMgcHJlZGljdCBpbmdyb3VwIGNob2ljZXMsIGNvbnRyb2xsaW5nIGZvciBzZWxmLWRlc2NyaXB0aXZlbmVzcz8KCmBgYHtyfQptIDwtIGdsbWVyKCBpbmdDaG9pY2VOIH4gc2NhbGUobmxzbG9wZSkgKyBzZWxmUmVzcC5aICsgICggc2NhbGUobmxzbG9wZSkgKyBzZWxmUmVzcC5aIHwgc3ViSUQpICsgKDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKdGlkeShtLGNvbmYuaW50PVRSVUUsZXhwb25lbnRpYXRlPVRSVUUsZWZmZWN0cz0iZml4ZWQiKQpgYGAKCiMgRG9lcyBhIGxpbmVhciB0cmVuZCBvZiBzaW1pbGFyaXR5LWJhc2VkIHByb2JhYmlsaXRpZXMgcHJlZGljdCBpbmdyb3VwIGNob2ljZXM/CgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IHNsb3BlLlogICsgcHJlZGljdGVkLlogKyAoIHNsb3BlLlogKyBwcmVkaWN0ZWQuWiB8IHN1YklEKSArICgxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCnRpZHkobSxjb25mLmludD1UUlVFLGV4cG9uZW50aWF0ZT1UUlVFLGVmZmVjdHM9ImZpeGVkIikKYGBgCgojIEJhY2t3YXJkcyBzb2x1dGlvbjogQ2FuIHlvdSBwcmVkaWN0IHNlbGYtZXZhbHVhdGlvbnMgZnJvbSBzaW1pbGFyaXR5IHRvIGluZ3JvdXAgYW5kIG91dGdyb3VwIGNob2ljZXM/CgpgYGB7cn0KbSA8LSBsbWVyKCBzY2FsZShzZWxmUmVzcCkgfiBzY2FsZShpbkdzaW0pICsgc2NhbGUob3V0R3NpbSkgKyAoICBzY2FsZShpbkdzaW0pICsgc2NhbGUob3V0R3NpbSkgfCBzdWJJRCkgKyAoMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUcmFpbikKc3VtbWFyeShtKQp0aWR5KG0sY29uZi5pbnQ9VFJVRSxlZmZlY3RzPSJmaXhlZCIpCmdncHJlZGljdChtLCBjKCJpbkdzaW0iKSkgJT4lIHBsb3Qoc2hvdy50aXRsZT1GKSArIHhsYWIoIlNpbWlsYXJpdHkgdG8gSW5ncm91cCBDaG9pY2VzIikgKyB5bGFiKCJTZWxmLUV2YWx1YXRpb24iKSArIGp0b29sczo6dGhlbWVfYXBhKCkKYGBgCgpgYGB7cn0KbSA8LSBnbG1lciggaW5nQ2hvaWNlTiB+IHNjYWxlKHNsb3BlKSAqIG5vdmVsICArICggc2NhbGUoc2xvcGUpICsgbm92ZWwgfCBzdWJJRCkgKyAoMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQp0aWR5KG0sY29uZi5pbnQ9VFJVRSxleHBvbmVudGlhdGU9VFJVRSxlZmZlY3RzPSJmaXhlZCIpCmBgYAoKYGBge3J9Cm0gPC0gZ2xtZXIoIGluZ0Nob2ljZU4gfiBlU0UgKyAoIGVTRSB8IHN1YklEKSAsIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKYGBgCgoKYGBge3J9Cm0gPC0gZ2xtZXIoIGluZ0Nob2ljZU4gfiBzU0UgKyAoIHNTRSB8IHN1YklEKSArICgxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCnRpZHkobSxjb25mLmludD1UUlVFLGV4cG9uZW50aWF0ZT1UUlVFLGVmZmVjdHM9ImZpeGVkIikKcjJiZXRhKG0pCmBgYAoKYGBge3J9Cm0gPC0gZ2xtZXIoIGluZ0Nob2ljZU4gfiBTRS5aICogbm92ZWwgKyAoIFNFLlogKyBub3ZlbCB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQp0aWR5KG0sY29uZi5pbnQ9VFJVRSxleHBvbmVudGlhdGU9VFJVRSxlZmZlY3RzPSJmaXhlZCIpCnIyYmV0YShtKQpgYGAKCgpgYGB7cn0KbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gU0UuWiArIHNjYWxlKGRlc2lyYWJpbGl0eSkgKyAoIFNFLlogfCBzdWJJRCkgKyAoIDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNjYWxlKG9TRSkgKyAoIHNjYWxlKG9TRSkgfCBzdWJJRCkgKyAoMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2NhbGUoaVNFKSArICggc2NhbGUoaVNFKSB8IHN1YklEKSArICgxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCgptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBTRS5aKmFzLmZhY3Rvcihub3ZlbCkgKyAoIFNFLlogKyBhcy5mYWN0b3Iobm92ZWwpIHwgc3ViSUQpICsgKCAgMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gU0UuWiphcy5mYWN0b3Iobm92ZWwpICsgc2NhbGUoZGVzaXJhYmlsaXR5KSArICggU0UuWiArIGFzLmZhY3Rvcihub3ZlbCkgKyBzY2FsZShkZXNpcmFiaWxpdHkpIHwgc3ViSUQpICsgKCAgMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiU0UiLCAibm92ZWwiKSkgJT4lIHBsb3QoKQpgYGAKCmBgYHtyfQptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBwcmVkaWN0ZWQuWipzY2FsZShSU0UpICsgKCBwcmVkaWN0ZWQuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygicHJlZGljdGVkLloiLCAiUlNFIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHByZWRpY3RlZC5aKnNjYWxlKFNDQykgKyAoIHByZWRpY3RlZC5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJwcmVkaWN0ZWQuWiIsICJTQ0MiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gcHJlZGljdGVkLloqc2NhbGUoRFMpICsgKCBwcmVkaWN0ZWQuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygicHJlZGljdGVkLloiLCAiRFMiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gcHJlZGljdGVkLloqc2NhbGUoTkZDKSArICggcHJlZGljdGVkLlogfCBzdWJJRCkgKyAoIDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoInByZWRpY3RlZC5aIiwgIk5GQyIpKSAlPiUgcGxvdCgpCgptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBwcmVkaWN0ZWQuWipzY2FsZShTSU5HLkluZCkgKyAoIHByZWRpY3RlZC5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJwcmVkaWN0ZWQuWiIsICJTSU5HLkluZCIpKSAlPiUgcGxvdCgpCgptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBwcmVkaWN0ZWQuWipzY2FsZShTSU5HLkludGVyKSArICggcHJlZGljdGVkLlogfCBzdWJJRCkgKyAoIDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoInByZWRpY3RlZC5aIiwgIlNJTkcuSW50ZXIiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gcHJlZGljdGVkLloqc2NhbGUoUHJvdG8pICsgKCBwcmVkaWN0ZWQuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygicHJlZGljdGVkLloiLCAiUHJvdG8iKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gcHJlZGljdGVkLloqc2NhbGUoU0kpICsgKCBwcmVkaWN0ZWQuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygicHJlZGljdGVkLloiLCAiU0kiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gcHJlZGljdGVkLloqc2NhbGUoTlRCKSArICggcHJlZGljdGVkLlogfCBzdWJJRCkgKyAoIHByZWRpY3RlZC5aIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJwcmVkaWN0ZWQuWiIsICJOVEIiKSkgJT4lIHBsb3QoKQpgYGAKCmBgYHtyfQptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBlbnRyb3B5Lloqc2NhbGUoUlNFKSArICggZW50cm9weS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJlbnRyb3B5LloiLCAiUlNFIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IGVudHJvcHkuWipzY2FsZShTQ0MpICsgKCBlbnRyb3B5LlogfCBzdWJJRCkgKyAoIDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoImVudHJvcHkuWiIsICJTQ0MiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gZW50cm9weS5aKnNjYWxlKERTKSArICggZW50cm9weS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJlbnRyb3B5LloiLCAiRFMiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gZW50cm9weS5aKnNjYWxlKE5GQykgKyAoIGVudHJvcHkuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiZW50cm9weS5aIiwgIk5GQyIpKSAlPiUgcGxvdCgpCgptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBlbnRyb3B5Lloqc2NhbGUoU0lORy5JbmQpICsgKCBlbnRyb3B5LlogfCBzdWJJRCkgKyAoIDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoImVudHJvcHkuWiIsICJTSU5HLkluZCIpKSAlPiUgcGxvdCgpCgptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBlbnRyb3B5Lloqc2NhbGUoU0lORy5JbnRlcikgKyAoIGVudHJvcHkuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiZW50cm9weS5aIiwgIlNJTkcuSW50ZXIiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gZW50cm9weS5aKnNjYWxlKFByb3RvKSArICggZW50cm9weS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJlbnRyb3B5LloiLCAiUHJvdG8iKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gZW50cm9weS5aKnNjYWxlKFNJKSArICggZW50cm9weS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJlbnRyb3B5LloiLCAiU0kiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gZW50cm9weS5aKnNjYWxlKE5UQikgKyAoIGVudHJvcHkuWiB8IHN1YklEKSArICggZW50cm9weS5aIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJlbnRyb3B5LloiLCAiTlRCIikpICU+JSBwbG90KCkKYGBgCgpgYGB7cn0KbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2xvcGUuWipzY2FsZShSU0UpICsgKCBzbG9wZS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJzbG9wZS5aIiwgIlJTRSIpKSAlPiUgcGxvdCgpCgptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBzbG9wZS5aKnNjYWxlKFNDQykgKyAoIHNsb3BlLlogfCBzdWJJRCkgKyAoIDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoInNsb3BlLloiLCAiU0NDIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNsb3BlLloqc2NhbGUoRFMpICsgKCBzbG9wZS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJzbG9wZS5aIiwgIkRTIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNsb3BlLloqc2NhbGUoTkZDKSArICggc2xvcGUuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygic2xvcGUuWiIsICJORkMiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2xvcGUuWipzY2FsZShTSU5HLkluZCkgKyAoIHNsb3BlLlogfCBzdWJJRCkgKyAoIDEgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoInNsb3BlLloiLCAiU0lORy5JbmQiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2xvcGUuWipzY2FsZShTSU5HLkludGVyKSArICggc2xvcGUuWiB8IHN1YklEKSArICggMSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygic2xvcGUuWiIsICJTSU5HLkludGVyIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNsb3BlLloqc2NhbGUoUHJvdG8pICsgKCBzbG9wZS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJzbG9wZS5aIiwgIlByb3RvIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNsb3BlLloqc2NhbGUoU0kpICsgKCBzbG9wZS5aIHwgc3ViSUQpICsgKCAxIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJzbG9wZS5aIiwgIlNJIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNsb3BlLloqc2NhbGUoTlRCKSArICggc2xvcGUuWiB8IHN1YklEKSArICggc2xvcGUuWiB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygic2xvcGUuWiIsICJOVEIiKSkgJT4lIHBsb3QoKQpgYGAKCmBgYHtyfQptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBzY2FsZShkZXNpcmFiaWxpdHkpKnNjYWxlKFJTRSkgKyAoIHNjYWxlKGRlc2lyYWJpbGl0eSkgfCBzdWJJRCkgKyAoIHNjYWxlKGRlc2lyYWJpbGl0eSkgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoImRlc2lyYWJpbGl0eSIsICJSU0UiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2NhbGUoZGVzaXJhYmlsaXR5KSpzY2FsZShTQ0MpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgc3ViSUQpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJkZXNpcmFiaWxpdHkiLCAiU0NDIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNjYWxlKGRlc2lyYWJpbGl0eSkqc2NhbGUoRFMpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgc3ViSUQpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJkZXNpcmFiaWxpdHkiLCAiRFMiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2NhbGUoZGVzaXJhYmlsaXR5KSpzY2FsZShORkMpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgc3ViSUQpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJkZXNpcmFiaWxpdHkiLCAiTkZDIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNjYWxlKGRlc2lyYWJpbGl0eSkqc2NhbGUoU0lORy5JbmQpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgc3ViSUQpICsgKCBzY2FsZShkZXNpcmFiaWxpdHkpIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJkZXNpcmFiaWxpdHkiLCAiU0lORy5JbmQiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2NhbGUoZGVzaXJhYmlsaXR5KSpzY2FsZShTSU5HLkludGVyKSArICggc2NhbGUoZGVzaXJhYmlsaXR5KSB8IHN1YklEKSArICggc2NhbGUoZGVzaXJhYmlsaXR5KSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiZGVzaXJhYmlsaXR5IiwgIlNJTkcuSW50ZXIiKSkgJT4lIHBsb3QoKQoKbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gc2NhbGUoZGVzaXJhYmlsaXR5KSpzY2FsZShQcm90bykgKyAoIHNjYWxlKGRlc2lyYWJpbGl0eSkgfCBzdWJJRCkgKyAoIHNjYWxlKGRlc2lyYWJpbGl0eSkgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoImRlc2lyYWJpbGl0eSIsICJQcm90byIpKSAlPiUgcGxvdCgpCgptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBzY2FsZShkZXNpcmFiaWxpdHkpKnNjYWxlKFNJKSArICggc2NhbGUoZGVzaXJhYmlsaXR5KSB8IHN1YklEKSArICggc2NhbGUoZGVzaXJhYmlsaXR5KSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiZGVzaXJhYmlsaXR5IiwgIlNJIikpICU+JSBwbG90KCkKCm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHNjYWxlKGRlc2lyYWJpbGl0eSkqc2NhbGUoTlRCKSArICggc2NhbGUoZGVzaXJhYmlsaXR5KSB8IHN1YklEKSArICggc2NhbGUoZGVzaXJhYmlsaXR5KSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiZGVzaXJhYmlsaXR5IiwgIk5UQiIpKSAlPiUgcGxvdCgpCmBgYAoKYGBge3J9Cm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IHByZWRpY3RlZC5aKm5vdmVsKnNjYWxlKE5GQykgKyAoIHByZWRpY3RlZC5aK25vdmVsIHwgc3ViSUQpICsgKCBTRS5aKmFzLmZhY3Rvcihub3ZlbCkgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoIlNFIiwgIm5vdmVsIiwiUlNFIikpICU+JSBwbG90KCkKYGBgCgpgYGB7cn0KbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gU0UuWiphcy5mYWN0b3Iobm92ZWwpKnNjYWxlKFNDQykgKyAoIFNFLlorYXMuZmFjdG9yKG5vdmVsKSB8IHN1YklEKSArICggU0UuWithcy5mYWN0b3Iobm92ZWwpIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJTRSIsIm5vdmVsIiwiU0NDIikpICU+JSBwbG90KCkKYGBgCgpgYGB7cn0KbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gU0UuWiphcy5mYWN0b3Iobm92ZWwpKnNjYWxlKERTKSArICggU0UuWithcy5mYWN0b3Iobm92ZWwpIHwgc3ViSUQpICsgKCBTRS5aK2FzLmZhY3Rvcihub3ZlbCkgfCB0cmFpdCksIGRhdGEgPSBmdWxsVGVzdCwgZmFtaWx5ID0gYmlub21pYWwsIGNvbnRyb2wgPSBnbG1lckNvbnRyb2wob3B0aW1pemVyID0gImJvYnlxYSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdEN0cmwgPSBsaXN0KG1heGZ1biA9IDEwMDAwMCkpLAogICAgbkFHUSA9IDEpCnN1bW1hcnkobSkKZ2dwcmVkaWN0KG0sIGMoIlNFIiwibm92ZWwiLCJEUyIpKSAlPiUgcGxvdCgpCmBgYAoKYGBge3J9Cm0gPC0gZ2xtZXIoIGFzLmZhY3RvcihpbmdDaG9pY2VOKSB+IFNFLloqYXMuZmFjdG9yKG5vdmVsKSpzY2FsZShORkMpICsgKCBTRS5aK2FzLmZhY3Rvcihub3ZlbCkgfCBzdWJJRCkgKyAoIFNFLlorYXMuZmFjdG9yKG5vdmVsKSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiU0UiLCJub3ZlbCIsICJORkMiKSkgJT4lIHBsb3QoKQpgYGAKCmBgYHtyfQptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBTRS5aKmFzLmZhY3Rvcihub3ZlbCkqc2NhbGUoU0lORy5JbmQpICsgKCBTRS5aK2FzLmZhY3Rvcihub3ZlbCkgfCBzdWJJRCkgKyAoIFNFLlorYXMuZmFjdG9yKG5vdmVsKSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiU0UiLCJub3ZlbCIsIlNJTkcuSW5kIikpICU+JSBwbG90KCkKYGBgCgpgYGB7cn0KbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gU0UuWiphcy5mYWN0b3Iobm92ZWwpKnNjYWxlKFNJTkcuSW50ZXIpICsgKCBTRS5aK2FzLmZhY3Rvcihub3ZlbCkgfCBzdWJJRCkgKyAoIFNFLlorYXMuZmFjdG9yKG5vdmVsKSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiU0UiLCJub3ZlbCIsIlNJTkcuSW50ZXIiKSkgJT4lIHBsb3QoKQpgYGAKCmBgYHtyfQptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBTRS5aKmFzLmZhY3Rvcihub3ZlbCkqc2NhbGUoUHJvdG8pICsgKCBTRS5aK2FzLmZhY3Rvcihub3ZlbCkgfCBzdWJJRCkgKyAoIFNFLlorYXMuZmFjdG9yKG5vdmVsKSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiU0UiLCAiUHJvdG8iKSkgJT4lIHBsb3QoKQpgYGAKCmBgYHtyfQptIDwtIGdsbWVyKCBhcy5mYWN0b3IoaW5nQ2hvaWNlTikgfiBTRS5aKmFzLmZhY3Rvcihub3ZlbCkqc2NhbGUoU0kpICsgKCBTRS5aK2FzLmZhY3Rvcihub3ZlbCkgfCBzdWJJRCkgKyAoIFNFLlorYXMuZmFjdG9yKG5vdmVsKSB8IHRyYWl0KSwgZGF0YSA9IGZ1bGxUZXN0LCBmYW1pbHkgPSBiaW5vbWlhbCwgY29udHJvbCA9IGdsbWVyQ29udHJvbChvcHRpbWl6ZXIgPSAiYm9ieXFhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3B0Q3RybCA9IGxpc3QobWF4ZnVuID0gMTAwMDAwKSksCiAgICBuQUdRID0gMSkKc3VtbWFyeShtKQpnZ3ByZWRpY3QobSwgYygiU0UiLCJub3ZlbCIsIlNJIikpICU+JSBwbG90KCkKYGBgCgpgYGB7cn0KbSA8LSBnbG1lciggYXMuZmFjdG9yKGluZ0Nob2ljZU4pIH4gU0UuWiphcy5mYWN0b3Iobm92ZWwpKnNjYWxlKE5UQikgKyAoIFNFLlorYXMuZmFjdG9yKG5vdmVsKSB8IHN1YklEKSArICggU0UuWithcy5mYWN0b3Iobm92ZWwpIHwgdHJhaXQpLCBkYXRhID0gZnVsbFRlc3QsIGZhbWlseSA9IGJpbm9taWFsLCBjb250cm9sID0gZ2xtZXJDb250cm9sKG9wdGltaXplciA9ICJib2J5cWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBvcHRDdHJsID0gbGlzdChtYXhmdW4gPSAxMDAwMDApKSwKICAgIG5BR1EgPSAxKQpzdW1tYXJ5KG0pCmdncHJlZGljdChtLCBjKCJTRSIsICJub3ZlbCIsICJOVEIiKSkgJT4lIHBsb3QoKQpgYGAKCgo=